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We present the most significant measurement of the cosmic microwave background (CMB) lensing potential to date (at a level of 40 cr), using 
temperature and polarization data from the Planck 2015 full-mission release. Using a polarization-only estimator, we detect lensing at a signifi¬ 
cance of 5 cr. We cross-check the accuracy of our measurement using the wide frequency coverage and complementarity of the temperature and 
polarization measurements. Public products based on this measurement include an estimate of the lensing potential over approximately 70% of 
the sky, an estimate of the lensing potential power spectrum in bandpowers for the multipole range 40<L<400, and an associated likelihood for 
cosmological parameter constraints. We find good agreement between our measurement of the lensing potential power spectrum and that found in 
the ACDM model that best fits the Planck temperature and polarization power spectra. Using the lensing likelihood alone we obtain a percent-level 
measurement of the parameter combination crgfl^ 25 = 0.591 ± 0.021. We combine our determination of the lensing potential with the £-mode 
polarization, also measured by Planck, to generate an estimate of the lensing B-mode. We show that this lensing B-mode estimate is correlated 
with the B-modes observed directly by Planck at the expected level and with a statistical significance of 10 cr, confirming Planck's sensitivity to 
this known sky signal. We also correlate our lensing potential estimate with the large-scale temperature anisotropies, detecting a cross-correlation 
at the 3 cr level, as expected because of dark energy in the concordance ACDM model. 


1. Introduction 

The cosmic microwave background (CMB) gives us a direct 
measurement of the early Universe when it first became trans¬ 
parent to radiation just 375 000 years after the Big Bang. It con¬ 
tains distinct signatures of the later Universe as well, which were 
imprinted by the process of gravitational lensing. The CMB pho- 
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tons that we observe today last scattered approximately 14 bil¬ 
lion years ago, and have travelled most of the way across the ob¬ 
servable Universe to reach us. During their journey, their paths 
were distorted by the gravitational tug of intervening matter, a 
subtle effect that may be measured statistically with high angu¬ 
lar resolution, low-noise observations of the CMB, such as those 
provided by Planck. In this paper, we present lensing measure¬ 
ments that are based on the full-mission 2015 data release. We 
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produce the most powerful measurement of CMB lensing to date 
with a 2.5 % constraint on the amplitude of the lensing potential 
power spectrum (or alternatively, a 40 cr detection of lensing ef¬ 
fects). 

The effect of lensing is to remap the CMB fluctuations, so 
that the observed anisotropy in direction h is in fact the unlensed, 
‘primordial’ anisotropy in the direction h + V0(/i), where <f>(h) 
is the CMB lensing potential defined by (e.g. Lewis & Challinor 
2006) 


<P(n) = -2 f dx Pa ~x)- (1) 

Jo Jk(x*)Jk(x) 

Here, x is conformal distance (with x* ~ 14000Mpc denot¬ 
ing the distance to the CMB last-scattering surface). The 
angular-diameter distance fx(x) depends on the curvature of the 
Universe, and is given by 


K~ l/2 sin (K 1/2 x) 


f K (x) = 


X 


|/cr 1/2 sinh(|/C| I/2 ^) 


for K > 0 (closed), 
for K - 0 (flat), 
for K < 0 (open). 


( 2 ) 


Finally, X ¥(xh; rf) is the (Weyl) gravitational potential at confor¬ 
mal distance x along the direction h at conformal time // (the 
conformal time today is denoted as 770 ). The lensing potential 
is an integrated measure of the mass distribution back to the 
last-scattering surface. The power spectrum of the lensing po¬ 
tential probes the matter power spectrum, which is sensitive to 
‘late-time’ parameters that modify the growth of structure such 
as neutrino mass (Smith et al. 2009). The amplitude of lensing 
effects is also a sensitive probe of geometrical parameters, such 
as the curvature of the Universe. 

The lens-induced remapping imprints distinctive statistical 
signatures onto the observed CMB fluctuations, which can be 
mined for cosmological information in a process known as 
lens reconstruction (Okamoto & Hu 2003). The past several 
years have been seen dramatic improvements in CMB lensing 
measurements, moving from first detections in cross-correlation 
(Smith et al. 2007; Hirata et al. 2008) to cosmologically use¬ 
ful measurements of the lensing potential power spectrum (Das 
et al. 2011; van Engelen et al. 2012; Planck Collaboration XVII 
2014). Recently, ground-based experiments have been able to 
detect the effects of lensing in polarization data as well (Hanson 
et al. 2013; POLARBEAR Collaboration 2014a,b; van Engelen 
et al. 2015). 

The results in this paper extend our earlier results in 
Planck Collaboration XVII (2014). These were based on the 
March 2013 Planck 1 nominal-mission data release, which con¬ 
tains approximately 15 months of temperature data alone. The 
additional information in the full-mission dataset, with 30 
months of Planck HFI temperature and polarization data (Planck 
Collaboration VII 2016; Planck Collaboration VIII 2016), al¬ 
lows us to improve our reconstruction noise levels by roughly 
a factor of two. Approximately half of this improvement comes 
from lower noise levels in temperature, and the other half from 
inclusion of polarization data. The improved lensing map is in¬ 
cluded as part of the Planck 2015 public data release (Planck 

* Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci¬ 
entific consortia funded by ESA member states and led by Principal 
Investigators from France and Italy, telescope reflectors provided 
through a collaboration between ESA and a scientific consortium led 
and funded by Denmark, and additional contributions from NASA 
(USA). 


Collaboration I 2016), as well as an estimate of the lensing po¬ 
tential power spectrum and associated likelihoods. In this paper 
we describe the creation of these products, as well as first sci¬ 
ence results based on them. We highlight the following science 
results. 

• We detect lensing /(-modes in the Planck data at a signifi¬ 
cance of 10 cr, using both a cross-correlation approach with 
the cosmic infrared background (CIB) as a tracer of the lens¬ 
ing potential, as well as a CMB-only approach using the 
TTEB trispectrum. This provides an important confirmation 
that Planck is sensitive to this known source of /(-modes on 
intermediate and small scales. 

• We make an improved measurement of the 3-point function 
(bispectrum) that is induced in the CMB by the correlation 
between lensing and the integrated Sachs-Wolfe (ISW) ef¬ 
fect, the latter being sourced by late-time acceleration. The 
lensing-ISW bispectrum is now detected at the 3 cr level. 

• Using only lensing information (along with well-motivated 
priors), we constrain the parameter combination crgQfJ ; 25 to 
roughly 3 %. 

The organization of this paper is as follows. In Sect. 2 we 
give a summary of our analysis pipeline for the full-mission 
data. Our analysis is very similar to that presented in Planck 
Collaboration XVII (2014), with straightforward extensions to 
polarization data, as well as a few additional improvements, 
and so we have deferred most technical details to appendices. 
In Sect. 3 we present our main results, including the lensing 
potential map, bandpower estimates of the lensing potential 
power spectrum, and implications for cosmological parameters. 
In Sect. 4 we present a suite of consistency and null tests to ver¬ 
ify our bandpower estimates, and in Sect. 5 we conclude. A se¬ 
ries of appendices provide more technical details of our analysis 
pipeline, modeling of the CIB, and dependence of the lensing 
power spectrum on ACDM parameters. 

2. Data and methodology 

Here we give a brief overview of the procedure that we use to 
measure the CMB lensing potential and its power spectrum from 
the Planck maps. We defer the detailed technical aspects of our 
pipeline to Appendix A. 

Our main results are based on a foreground-cleaned map of 
the CMB synthesized from the raw Planck 2015 full-mission 
frequency maps using the SMICA code (Planck Collaboration 
IX 2016). This foreground-cleaned map combines all nine fre¬ 
quency bands from 30 GHz to 857 GHz with scale-dependent 
coefficients chosen to provide unit response to the CMB with 
minimal variance; for more details see Planck Collaboration IX 
(2016)*. On the small scales where lensing effects are most im¬ 
portant (multipoles above i = 1000), most of the CMB informa¬ 
tion in the foreground-cleaned maps originates with the 143 GHz 
and 217 GHz channel data. These channels have beams that 
are approximately Gaussian with full-width-at-half-maximum 
(FWHM) parameters of 7' and 5', respectively. Their noise is 
approximately white for 1000 < l < 3000, with levels of 
30 fjK arcrnin (60 fjK. arc min) in temperature (polarization) at 
143 GHz, and 40 /rK arcrnin (95 /jK arcmin) at 217 GHz. 

We reconstruct the CMB lensing potential using quadratic 
estimators that exploit the statistical anisotropy induced by lens¬ 
ing, following Okamoto & Hu (2003). Neglecting the lensing of 

• We use a slightly earlier version of SMICA maps than those in the 
public release. They differ by pixel-smoothing and missing-pixel mask, 
however this has no significant effect on the results. 
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primordial B-modes there are five possible estimators, denoted 
by cj) TT , cj) TE , (p EE , 4 > eb , and <p TB , which are based on various cor¬ 
relations of the CMB temperature ( T ) and polarization (E and 
B ). In addition, we can form a minimum-variance estimator that 
combines all five estimators, which we denote as $> MV . In Fig. 1 
we plot the lens reconstruction noise levels for these estimators. 
The most powerful estimator is TT, although the TE and EE 
estimators are also useful on large angular scales. 



Fig. 1 Lens reconstruction noise levels /V‘ ‘ for the TT, TE, EE, 
EB, and TB estimators applied to the SMICA full-mission CMB 
map. The noise level for their minimum-variance combination 
(MV) is also shown. The fiducial ACDM theory power spectrum 
Cf lld used in our Monte Carlo simulations is plotted as the black 
solid line. 

The quadratic lensing estimators take inverse-variance fil¬ 
tered CMB multipoles as input. We obtain these using a filter that 
masks the Galaxy and point sources, and also bandpass filters 
the data in harmonic space to 100 < ( < 2048. For our baseline 
analysis, we start with the masks used in the analysis of Planck 
Collaboration XVII (2014), with a Galaxy mask that removes 

30.2 % of the sky and a point-source mask that removes an ad¬ 
ditional 0.7 % of the sky. This removes the brightest Sunyaev- 
Zel’dovich (SZ) clusters as described in Planck Collaboration 
XVII (2014); contamination of the lensing reconstruction from 
residual unresolved SZ is expected to be at the percent level 
for Planck (van Engelen et al. 2014), even without compo¬ 
nent separation. ’ This residual signal, and that from unresolved 
point-sources, is further reduced by a correction that we make 
to the power spectrum of the reconstructed lensing potential 
(see Appendix A. 3). Finally, we apply SMICA-specific temper¬ 
ature and polarization masks described in Planck Collaboration 
IX (2016). Combining all three sets of masks leaves a total of 

67.3 % of the sky for analysis. 


t We mask clusters with S /N > 5 in the Planck cluster catalogue 
that accompanied the 2013 release (PSZ1; Planck Collaboration XXIX 
2014). The mass limit is redshift dependent, with a high-redshift limit 
at 80% completeness of M 500 * 6 x 10 14 M G . If we conservatively 
adopt this mass limit (we certainly mask low-redshift clusters to lower 
masses), the results of van Engelen et al. (2014) suggest that the SZ 
trispectrum at 143 GHz gives a positive bias of at most a few percent 
in the reconstructed lensing power spectrum for multipoles L < 2000, 
while the lensing-SZ-SZ bispectrum gives a negative bias of similar 
size. The SZ signal is much smaller in the 217 GHz temperature data, 
which is around the null of the (thermal) SZ effect. 


We estimate the power spectrum of the lensing potential Cf' 
using the auto- and cross-spectra of the quadratic lensing esti¬ 
mators. 5 These spectra probe the 4-point function of the lensed 
CMB, specifically the connected (trispectrum) part of the 4-point 
function that is sourced by lensing. They also contain contribu¬ 
tions from the disconnected part of the 4-point function (which 
is non-zero even in the absence of lensing effects). We estimate 
this contribution and subtract it, as well as several other smaller 
bias terms, obtaining an estimate C^' of the power spectrum. 
This procedure is discussed in Appendix A. 3. For cosmologi¬ 
cal parameter constraints, we use a Gaussian log-likelihood in 
bandpowers of the estimated lensing power spectrum, given by 

- 2 log - Cf lh ) [x 1 1'' S^'(Cf? - cff’ th ). (3) 

Here, bins are indexed by i and j; is the bandpower bin¬ 
ning function for the ith bin, and X is a covariance matrix for 
the bin estimates. Paired upper/lower indices are summed over. 
The Cf th is the theoretical expectation value of the estimated 
Cf for the set of cosmological and nuisance parameters under 
consideration. This generally differs from the theory spectrum 
Cf* at the same cosmological parameters due to the way that our 
power spectrum estimates are normalized, and corrected for ad¬ 
ditional trispectrum couplings, with a fiducial model. Both the 
binning function and cf^’ th are discussed in Appendix C. The 
binning function is chosen to have unit response to a fiducial 
theory spectrum Cf f ’ A ' il , and so we denote 

if = S E Cf (4) 

as the amplitude of the power spectrum for a particular bin rela¬ 
tive to the fiducial expectation (with A = 1 for Cf^ = cf^’ fid ). 

To characterize the variance of our lensing potential esti¬ 
mates, as well as to estimate several bias terms, we use simulated 
Planck maps. These are based on the Full Focal Plane 8 (FFP8) 
Monte Carlo simulation set described in Planck Collaboration 
XII (2016). As discussed there, the Planck maps were effectively 
renormalized by approximately 2-3% in power in the time be¬ 
tween the generation of FFP8 and the final Planck full-mission 
maps. To account for this, we rescale the CMB component of the 
simulations by a factor of 1.0134 before analysis. The FFP8 sim¬ 
ulations do not include contributions from residual foregrounds 
(Galactic dust, as well as unmasked point sources), and also un¬ 
derestimate the noise power spectra by several percent at high-/ - . 
We account for this missing power simply by adding coloured 
Gaussian noise to the simulations to make their TT, EE, and 
BB power spectra agree with the data. This approach implicitly 
assumes that any non-Gaussianity of these residual components 
does not couple significantly to our lensing estimates. We per¬ 
form consistency tests in Sect. 4 to check the validity of these 
assumptions. 

Throughout this paper we use a spatially-flat fiducial cos¬ 
mology with baryon density given by a>\, = El^h 1 = 0.0222, 
cold dark matter density io c = El Q h 1 = 0.1203, neutrino en¬ 
ergy density io v = Q. v h 2 = 0.00064 (corresponding to two mass¬ 
less neutrinos and one massive with mass 0.06 eV), Hubble con¬ 
stant // (l = 100/ikms 1 Mpc 1 with h = 0.6712, spectral in¬ 
dex of the power spectrum of the primordial curvature perturba¬ 
tion n s = 0.96, amplitude of the primordial power spectrum (at 
k = 0.05 Mpc -1 ) A s = 2.09 x 10 -9 , and Thomson optical depth 

4 In this paper, we use multipole indices LM for the lens reconstruc¬ 
tion, reserving multipole indices £m for the CMB fields. 


3 












Planck Collaboration: Gravitational lensing by large-scale structures with Planck 


through reionization r = 0.065. These cosmological parame¬ 
ters also form the basis for the FFP8 Monte Carlo simulation 
set. In addition to rescaling the FFP8 maps as already discussed, 
we have also adjusted the power spectra of the fiducial model 
by rescaling the CMB temperature and polarization spectra by a 
factor of 1.0134 2 , and the temperature-lensing cross-correlation 
C ; / ' 2 by 1.0134. We have not applied any scaling to the fidu¬ 
cial lensing power spectmm. Our reconstruction methodology 
(in particular, the renormalization corrections, addition of fore¬ 
ground power, and realization-dependent bias corrections that 
we apply, discussed in Appendix C) renders the cosmological 
interpretation of our lensing estimates insensitive to errors in the 
fiducial model power spectra and simulations. 


3. Results 

In this section, we provide a summary of the first science re¬ 
sults obtained with the minimum-variance lens reconstruction 
from the Planck full-mission data. The lensing potential map is 
presented in Sect. 3.1, and this is combined in Sect. 3.2 with 
the E-mode polarization measured by Planck to obtain a map 
of the expected /l-mode polarization due to lensing. We fur¬ 
ther show that this is correlated with the /1-modes measured by 
Planck at the expected level. In Sect. 3.3, we cross-correlate the 
reconstructed lensing potential with the large-angle temperature 
anisotropies to measure the correlation sourced by the ISW 
effect. Finally, the power spectmm of the lensing potential is pre¬ 
sented in Sect. 3.4. We use the associated likelihood alone, and 
in combination with that constructed from the Planck temper¬ 
ature and polarization power spectra (Planck Collaboration XI 
2016), to constrain cosmological parameters in Sect. 3.5. 

3 . 1. Lensing potential 

In Fig. 2 we plot the Wiener-filtered minimum-variance lensing 
estimate, given by 


s~i(p(p, fid 

nWF _ C L 2MV 

VLM - 00 . fid 00 tLM ’ 


(5) 


where C ^’ fid is the lensing potential power spectrum in our fidu¬ 
cial model and A*' 4 is the noise power spectmm of the recon¬ 
struction. As we shall discuss in Sect. 4.5, the lensing potential 
estimate is unstable for L < 8, and so we have excluded those 
modes for all analyses in this paper, as well as in the MV lensing 
map. 

As a visual illustration of the signal-to-noise level in the lens¬ 
ing potential estimate, in Fig. 3 we plot a simulation of the MV 
reconstruction, as well as the input (f> realization used. The re¬ 
construction and input are clearly correlated, although the recon¬ 
struction has considerable additional power due to noise. As can 
be seen in Fig. 1, even the MV reconstruction only has S/N « 1 
for a few modes around L » 50. 

The MV lensing estimate in Fig. 2 forms the basis for a 
public lensing map that we provide to the community (Planck 
Collaboration I 2016). The raw lensing potential estimate has a 
very red power spectrum, with most of its power on large angular 
scales. This can cause leakage issues when cutting the map (for 
example to cross-correlate with an additional mass tracer over a 
small portion of the sky). The lensing convergence k defined by 


klm 


UL + 1 ) 

- 2 - Qlm, 


(6) 



0 WF (Data) 


Fig. 2 Lensing potential estimated from the SMICA full-mission 
CMB maps using the MV estimator. The power spectrum of 
this map forms the basis of our lensing likelihood. The estimate 
has been Wiener filtered following Eq. (5), and band-limited to 
8 < L < 2048. 



Fig. 3 Simulation of a Wiener-filtered MV lensing reconstruc¬ 
tion (upper) and the input </> realization (lower), filtered in the 
same way as the MV lensing estimate. The reconstruction and 
input are clearly correlated, although the reconstruction has con¬ 
siderable additional power due to noise. 


has a much whiter power spectmm, particularly on large angular 
scales. The reconstruction noise on k is approximately white as 
well (Bucher et al. 2012). For this reason, we provide a map 
of the estimated lensing convergence k rather than the lensing 
potential <p. 

3.2. Lensing B-mode power spectrum 

The odd-parity B-mode component of the CMB polarization is 
of great importance for early-universe cosmology. At first order 
in perturbation theory it is not sourced by the scalar fluctuations 
that dominate the temperature and polarization anisotropies, and 
so the observation of primordial /1-modes can be used as a 
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uniquely powerful probe of tensor (gravitational wave) or vec¬ 
tor perturbations in the early Universe. A detection of /(-mode 
fluctuations on degree angular scales, where the signal from 
gravitational waves is expected to peak, has recently been re¬ 
ported at 150 GHz by the BICEP2 collaboration (Ade et al. 

2014) . Following the joint analysis of BICEP2 and Keck Array 
data (also at 150 GHz) and the Planck polarization data, primar¬ 
ily at 353 GHz (BICEP2/Keck Array and Planck Collaborations 

2015) , it is now understood that the 3-mode signal detected 
by BICEP2 is dominated by Galactic dust emission. The joint 
analysis gives no statistically-significant evidence for primor¬ 
dial gravitational waves, and establishes a 95 % upper limit 
h ).05 < 0.12. This still represents an important milestone for 
/(-mode measurements, since the direct constraint from the /(- 
mode power spectrum is now as constraining as indirect, and 
model-dependent, constraints from the TT spectrum (Planck 
Collaboration XIII 2016). 

In addition to primordial sources, the effect of gravitational 
lensing also generates /(-mode polarization. The displacement of 
lensing mixes /(-mode polarization into /(-mode as (Smith et al. 
2009) 

B 5. - (-D~ Z Z (me -m B M ) ( 7 ) 

LM e E m E ' ' 

where is a weight function and the bracketed term is a 

Wigner-3 j symbol. On scales f'n < 1000 the lensing /(-mode 
power spectrum resembles that of white noise, with a level of 
about 5 pK arcmin. This lensing power acts as a potential source 
of confusion for the measurement of primordial /(-modes, which 
can be estimated and ultimately removed in a process of delens- 
ing. Given an estimate for the lensing potential (f> and the E-mode 
polarization measured by Planck we can synthesize a lensed B- 
mode map for this purpose using Eq. (7). The 5 fiK arcmin level 
of the lensing /(-mode power spectrum is an order of magnitude 
lower than the Planck 2015 noise levels, and so delensing does 
not significantly improve our /(-mode measurements; however, 
the cross-correlation of the lensing /(-mode template with the 
observed /(-mode sky provides a useful check on the ability of 
Planck to measure this known source of /(-modes. 

We show the results of such a cross-correlation in Fig. 4, 
finding good agreement with the expected lensing /(-mode 
power spectrum. In addition to our fiducial MV lensing potential 
estimate, which uses both temperature and polarization data, we 
have also estimated the lensing /(-mode power spectrum using 
the TT -only lensing estimator to measure </>, as well as the CIB 
fluctuations measured by the 545 GHz Planck channel . 1 We see 
good agreement in all cases with the expected power; constrain¬ 
ing the overall amplitude of the lensing /(-mode power spectrum 
A b (relative to the predicted spectrum in our fiducial model) for 
a large bin from 8 < /’b < 2048 we measure amplitudes of 

^2048 = 0.93 ± 0.08 (0 MV ), 

^ 8-72048 = 0.95 ± 0.09 ( cf > TT ), 

^8-^2048 = 0-93 ± 0.10 (CIB) 

for the three estimates, each corresponding to a roughly 10 cr de¬ 
tection of lensing /(-mode power in the Planck data. The shape of 
the cross-correlation is also in good agreement with expectation. 

t To calculate the scaling between the CIB map and the lensing po¬ 
tential (j>, we model the CIB using the simple model of Hall et al. (2010) 
for its redshift and frequency dependence. Further details are given in 
Appendix D. 


Taking the bins in Fig. 4 as independent, forming a x 1 relative 
to the theory model and comparing to the distribution from sim¬ 
ulations we obtain probability-to-exceed (PTE) values of 48%, 
71%, and 78% using the MV, TT, and CIB lensing estimates, 
respectively. 



Fig. 4 Measurement of the lensing /(-mode power spectrum, us¬ 
ing cross-correlation with estimates of the lensing potential as 
discussed in Sect. 3.2. The theoretical lensing /(-mode power 
spectrum, for the parameters of the fiducial cosmological model 
of Sect. 2, is plotted as dashed black. Green squares and blue 
triangles are results using the TT and MV </> reconstructions, re¬ 
spectively, to construct the lens-induced /(-mode template, while 
red circles use the CIB (from the Planck 545 GHz channel) to 
construct a proxy for 0. Lensing /(-mode power is detected with 
the expected scale dependence and amplitude at a significance 
level of approximately 10 cr. 


3.3. Lensing-ISW bispectrum 

As photons travel towards us from the last scattering surface, 
they are not only deflected by gravitational lensing, they also re¬ 
ceive net red/blueshifts from gravitational potentials that evolve 
if they are crossed at late times. This phenomenon, known as 
the integrated Sachs-Wolfe (ISW) effect, is believed to gen¬ 
erate anisotropies in the observed CMB temperature on large 
((' < 100) angular scales. It is of particular interest because the 
decay of gravitational potentials, which produces the ISW effect, 
does not occur during matter domination, but only at redshifts 
Z < 2 when dark energy becomes dynamically important. The 
ISW effect can be detected statistically by cross-correlating the 
observed temperature anisotropies with a tracer of structure at 
these redshifts. Here we use the lensing potential, which is well- 
matched to the ISW effect (Hu 2002). Current Planck results 
incorporating additional external tracers of large-scale structure 
are summarized in Planck Collaboration XXI (2016). 

In Fig. 5 we plot the cross-correlation between the 
MV lens reconstruction and the CMB temperature. This cross¬ 
correlation probes the bispectrum, or three-point correlation 
function of the CMB, which is due to the correlation of the lens¬ 
ing and ISW effects. The measurement is noisy, due to a com¬ 
bination of noise in the lens reconstruction and cosmic variance 
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in the temperature (which ultimately limits the detection of the 
lensing-ISW bispectrum to about 9 cr; Lewis et al. 201 1). 

To determine the overall detection significance for the cross¬ 
correlation, we use the minimum-variance bispectrum estimator 

, fTtp, fid - rj-<* 

XT<P _ 1 _Z,_ iLM. _ 1 LM /O', 

N T *{^ /sky (Cf M +Nt*)C T L T ’ M ' ’ 

where N T<fl is a normalization determined from simulations . 1 For 
the MV lens reconstruction, using 8 < L < 100 we measure an 
amplitude 

A ^ 100 = 0.90 + 0.28 (MV), (9) 

which is consistent with the theoretical expectation of unity and 
non-zero at just over 3 cr. Using the TT -only lensing estimate 
rather than the MV lensing estimate in the cross-correlation, we 
obtain 

A^i 00 = 0.68 ± 0.32 (TT). (10) 

Using simulations, we measure an rms difference between the 
TT and MV bispectrum amplitudes of 0.18 (roughly equal 
to the quadrature difference of their error bars, which is 
V0.32 2 - 0.28 2 = 0.15). Therefore the difference of amplitudes, 
AAg^oo = 0 . 22 , is compatible with the expected scatter. 



Fig. 5 Lensing-ISW bispectrum on large angular scales. The 
cross-spectrum between the MV lensing potential estimate and 
the temperature anisotropy is plotted for bins of width AL =15, 
covering the multipole range L = 8-98. The dashed line shows 
the predicted cross-spectrum in the fiducial model. The lensing- 
ISW bispectrum is detected at just over 3 cr significance. 


In Planck Collaboration XVII (2014), using the TT lensing 
estimator on the multipole range 10 < L < 100 we measured a 
somewhat higher value for the lensing-ISW bispectrum ampli¬ 
tude of A^ = 0.85 ± 0.35." We expect the difference with 

t We find N 1(4 is within 4% of the analytical expectation 



1 1 




* The amplitude quoted in Planck Collaboration XVII (2014) is ac¬ 
tually A^Uioc) 3 = 0-28 ± 0.32, however this is measured with respect 


respect to the 2015 TT measurement to have standard deviation 
of approximately Vo.35 2 - 0.32 2 = 0.14, and so the observed 
difference of 0.85 - 0.68 = 0.17 is reasonable. 

3.4. Lensing potential power spectrum 

In Fig. 6 we plot our estimate of the lensing potential power 
spectrum obtained from the MV reconstruction, as well as sev¬ 
eral earlier measurements. We see good agreement with the 
shape in the fiducial model, as well as earlier measurements (de¬ 
tailed comparisons with the 2013 spectrum are given later in this 
section). In Sect. 4 we perform a suite of internal consistency 
and null tests to check the robustness of our lensing spectrum to 
different analysis and data choices. 

We estimate the lensing potential power spectrum in band- 
powers for two sets of bins: a ‘conservative’ set of eight 
uniformly-spaced bins with AL = 45 in the range 40 < L < 400; 
and an ‘aggressive’ set of 18 bins that are uniformly spaced in 
L 06 over the multipole range 8 < L < 2048. The conservative 
bins cover a multipole range where the estimator signal-to-noise 
is greatest. They were used for the Planck 2013 lensing likeli¬ 
hood described in Planck Collaboration XVII (2014). The ag¬ 
gressive bins provide good sensitivity to the shape of the lensing 
power on both large and small scales, however they are more 
easily biased by errors in the mean-field corrections (which are 
large at L < 40) and the disconnected noise bias corrections (at 
L > 400). Results for the bandpower amplitudes Af, defined in 
Eq. (4), are given in Table 1 for both sets of multipole bins. 

Nearly all of the internal consistency tests that we present 
in Sect. 4 are passed at an acceptable level. There is, however, 
mild evidence for a correlated feature in the curl-mode null-test, 
centred around L ~ 500. The L range covered by this feature 
includes 638 < L < 762, for which the (gradient) lensing re¬ 
construction bandpower is 3.6 cr low compared to the predicted 
power of the fiducial model. In tests of the sensitivity of pa¬ 
rameter constraints to lensing multipole range, described later 
in Sect. 3.5.4, we find shifts in some parameters of around 1 cr 
in going from the conservative to aggressive range, with negli¬ 
gible improvement in the parameter uncertainty. Around half of 
these shifts come from the outlier noted above. For this reason, 
we adopt the conservative multipole range as our baseline here, 
and in other Planck 2015 papers, when quoting constraints on 
cosmological parameters. However, where we quote constraints 
on amplitude parameters in this paper, we generally give these 
for both the aggressive and conservative binning. The aggres¬ 
sive bins are also used for all of the C l! j l!l bandpower plots in this 
paper. 

Estimating an overall lensing amplitude (following Eq. 4) 
relative to our fiducial theoretical model, for a single bin over 
both the aggressive and conservative multipole ranges we find 

A^Xoo = 0 987 ± 0-025, (11a) 

A^ g = 0.983 ± 0.025. (lib) 

These measurements of the amplitude of the lensing power spec¬ 
trum both have precision of 2.5%, and are non-zero at 40 cr. 
Given the measured amplitude, it is clear that our overall lens¬ 
ing amplitude estimate is consistent with our fiducial ACDM 
model (which has A = 1). The shape of our measurement is also 
in reasonable agreement with the fiducial model. Marginalizing 

to a slightly different fiducial cosmology than the one used here. Those 
measurements have been renormalized to the fiducial model used for 
this paper with a factor of 1.09. 
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over an overall amplitude parameter, for the 19 bins plotted 
in Fig. 6 we obtain a x 2 with respect to our fiducial model 
of 28 (with AW = 18), with a corresponding PTE of 6%. A 
large portion of this x 2 is driven by the outlier bandpower with 
638 < L < 762. Removing this bandpower gives a x~ of 16.5 
(JVdof = 17), with a PTE of 49%. For the more conservative 
multipole range 40 < L < 400, using eight linear bins we obtain 
* 2 of 8.9 (JVdof = 7), with a PTE of 26 %. 

Finally, we note that the lensing bandpowers measured here 
are in good agreement with the Planck 2013 results (Planck 
Collaboration XVII 2014). The only clear visual difference be¬ 
tween the two measurements in Fig. 6 is for the L — 40- 
65 bin, where the 2013 bandpower is significantly higher than 
the current measurement. As previously discussed in Planck 
Collaboration XVII (2014), this bin amplitude was sensitive to 
the choice of foreground cleaning, and decreased when using 
a component-separated map for the reconstruction (as is done 
here), and so this difference is expected. To compare measure¬ 
ments in further detail, we calculate a x 2 f° r the difference be¬ 
tween our 2013 and 2015 MV lensing bandpower estimates as¬ 
suming a diagonal covariance as 


T2013-2015 


z 




Lb 12013 


c 


L b 12015 


r 


v “ r KU- v “ r K , LJ 


( 12 ) 


where Lb indexes the bandpower bins. We approximate the vari¬ 
ance of the differences of the bandpowers as the difference of 
the variances in forming this^f 2 . Using the conservative bins for 
40 < L < 400, which were used for the 2013 lensing likelihood, 
we obtain a value of T 2013-2015 = 10-8 (AW = 8), with a corre¬ 
sponding PTE of 22%. The 2013 result is temperature only; if 
we compare it directly to the 2015 TT lensing bandpowers, we 
obtain a value of x 2 = 7.2 (AW = 8, PTE = 52%). We note 
that some care should be taken when comparing the measured 
lensing amplitudes between 2013 and 2015, because the fidu¬ 
cial spectra against which they are measured differ significantly. 
In 2013, for example, we measured A^^ = 0.943 ± 0.040. 
However, this was measured with respect to a fiducial C' 1 " 1 ' that 
is between 6 % and 7 % higher than the one used here in this 
multipole range. Renormalizing the 2013 measurements to de¬ 
termine the amplitude for the 2015 fiducial cosmology, we ob¬ 
tain 1.005 + 0.043, which can be compared directly to the mea¬ 
surement of A^q^Jqq = 0.987 ± 0.025 for 2015. It is interesting 
to note that there is a 2 % difference in power between the Planck 
2013 and 2015 maps from HFI due to a change in the calibra¬ 
tion, as well as an improved beam model (Planck Collaboration 
VII 2016). This shift does not couple to the lensing reconstruc¬ 
tion; lensing is a geometric effect, and so an overall change in 
the calibration does not affect the lens reconstruction procedure, 
provided that a consistent power spectrum is used to analyse the 
maps, as is done here. 


3.5. Likelihoods and cosmological parameter constraints 

In this section we discuss cosmological implications of the 
Planck lensing potential power spectrum estimate. We use a 
Gaussian likelihood in the measured bandpowers, with correc¬ 
tions to account for errors in the fiducial CMB power spectra 
that are used to normalize and debias the lensing estimates. This 
procedure is described in more detail in Appendix C. 

In the following subsections, we discuss the information that 
can be gleaned from the lensing potential power spectrum alone, 
or in conjunction with the Planck 2015 CMB power spectrum 


Table 1 Lensing potential power spectrum estimates from the 
MV lens reconstruction. The A' 1 ’ values are dimensionless, with 
A <!> - 1 for a measured spectrum equal to the lensing potential 
power spectrum of our fiducial cosmology (described in Sect. 2). 
The final column lists bandpower estimates, averaged within 
each bin. 


Lensing power spectrum bandpowers 

^-Tnin 

^-Tnax 

10 7 [L(L+l)] 2 Cf/2 tt 

Conservative multipole range (40<L<400) 


40 

84 

1.01 ±0.05 

1.28 ±0.06 

85 

129 

1.02 ±0.04 

1.00 ±0.04 

130 

174 

0.93 ± 0.05 

0.70 ± 0.04 

175 

219 

0.89 ± 0.07 

0.52 ± 0.04 

220 

264 

0.86 ±0.10 

0.41 ± 0.05 

265 

309 

1.03 ±0.11 

0.41 ± 0.05 

310 

354 

1.24 ±0.13 

0.41 ± 0.04 

355 

400 

0.97 ±0.16 

0.27 ± 0.05 

Aggressive multipole range (8<L<2048) 


8 

20 

1.09 ±0.17 

1.29 ±0.21 

21 

39 

1.01 ±0.09 

1.39 ±0.13 

40 

65 

1.05 ±0.07 

1.40 ±0.09 

66 

100 

1.00 ±0.05 

1.13 ±0.06 

101 

144 

0.99 ± 0.04 

0.88 ± 0.04 

145 

198 

0.89 ± 0.05 

0.60 ± 0.04 

199 

263 

0.90 ± 0.08 

0.45 ± 0.04 

264 

338 

1.08 ±0.10 

0.40 ± 0.04 

339 

425 

1.08 ±0.11 

0.30 ± 0.03 

426 

525 

0.98 ±0.13 

0.20 ± 0.03 

526 

637 

0.69 ±0.17 

0.11 ±0.03 

638 

762 

0.18 ±0.23 

0.02 ± 0.03 

763 

901 

0.73 ± 0.27 

0.07 ± 0.02 

902 

1054 

1.11 ±0.40 

0.08 ± 0.03 

1055 

1221 

0.27 ± 0.55 

0.01 ± 0.03 

1222 

1404 

1.06 ±0.80 

0.05 ± 0.03 

1405 

1602 

0.91 ± 1.18 

0.03 ± 0.04 

1603 

1816 

-2.03 ± 1.75 

-0.06 ± 0.05 

1817 

2048 

-0.54 ±2.16 

-0.01 ± 0.05 


likelihood. Some of the results in this section summarize the 
parts of Planck Collaboration XIII (2016) that are directly re¬ 
lated to lensing. 


3.5.1. Constraints from CMB lensing alone 

The measurement of the lensing power spectrum is obtained 
from the 4-point function of the observed CMB anisotropies, 
and hence depends in general on both the lensing and CMB 
power spectra (see Appendix C). However the latter are now 
well measured by Planck , and there is only a weak dependence 
on cosmological parameters given the observed spectra (mainly 
via different cosmological models changing foreground param¬ 
eters, and hence the amplitude of underlying CMB power spec¬ 
tra). In this section we fix the CMB power spectra to a ACDM 
Planck TT+lowP best fit (i.e. using TT on all scales and low-f’ 
polarization), and consider the conditional probability of differ¬ 
ent cosmological parameters given only the lensing reconstruc¬ 
tion power spectrum. 

To avoid marginalizing over very unrealistic values of poorly 
constrained parameters, we adopt several well-motivated priors 
when considering constraints from lensing alone. 
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Fig. 6 Planck 2015 full-mission MV lensing potential power spectrum measurement, as well as earlier measurements using the 
Planck 2013 nominal-mission temperature data (Planck Collaboration XVII 2014), the South Pole Telescope (SPT, van Engelen 
et al. 2012), and the Atacama Cosmology Telescope (ACT, Das et al. 2014). The fiducial ACDM theory power spectrum based on 
the parameters given in Sect. 2 is plotted as the black solid line. 


• The optical depth to reionization is fixed to r = 0.07, because 
lensing deflections are independent of reionization (and scat¬ 
tering and subsequent lensing from sources at reionization is 
negligible). 

• The baryon density is given a Gaussian 1 cr prior Q h /z 2 = 
0.0223 ± 0.0009, as measured independently from big bang 
nucleosynthesis models combined with quasar absorption 
line observations (Pettini & Cooke 2012). 

• The scalar spectral index is given a broad prior n s = 0.96 ± 
0.02; results are only weakly sensitive to this choice, within 
plausible bounds. 

• A top-hat prior is used for the reduced Hubble constant, 
0.4 < lr < 1. This limits the extent of the parameter degener¬ 
acy, but does not affect the results over the region of interest 
for joint constraints. 

In addition to the priors above, we adopt the same sampling 
priors and methodology as Planck Collaboration XIII (2016), 1 
using CosmoMC and camb for sampling and theoretical predic¬ 
tions (Lewis & Bridle 2002; Lewis et al. 2000). In the ACDM 
model, as well as Cl\,h 2 and n s , we sample A s , O c /z 2 , and the 
(approximate) acoustic-scale parameter $mc- Alternatively, we 
can think of our lensing-only results as constraining the sub¬ 
space of Q m , //(,, and cr s . Figure 7 shows the corresponding 
constraints from CMB lensing, along with tighter constraints 
from combining with additional external baryon acoustic oscil¬ 
lation (BAO) data, compared to the constraints from the Planck 
CMB power spectra. The contours overlap in a region of accept¬ 
able Hubble constant values, and hence are compatible. To show 

t For example, we split the neutrino component into approximately 
two massless neutrinos and one with 2 m v = 0.06 eV, by default. 


the multi-dimensional overlap region more clearly, the red con¬ 
tours show the lensing constraint when restricted to a reduced- 
dimensionality space with 0 mc fixed to the value accurately mea¬ 
sured by the CMB power spectra; the intersection of the red and 
black contours gives a clearer visual indication of the consis¬ 
tency region in the O m -rr 8 plane. 

The lensing-only constraint defines a band in the D m -cr 8 
plane, with the well-constrained direction corresponding ap¬ 
proximately to the constraint 

ct 8 Q^ 25 = 0.591 +0.021 (lensing only; 68 %). (13) 

This parameter combination is measured with approximately 
3.5% precision. 

The dependence of the lensing potential power spectrum on 
the parameters of the ACDM model is discussed in detail in 
Appendix E; see also Pan et al. (2014). Here, we aim to use 
simple physical arguments to understand the parameter degen¬ 
eracies of the lensing-only constraints. In the flat ACDM model, 
the bulk of the lensing signal comes from high redshift (z > 0.5) 
where the Universe is mostly matter-dominated (so potentials are 
nearly constant), and from lenses that are still nearly linear. For 
fixed CMB (monopole) temperature, baryon density, and n s , in 
the ACDM model the broad shape of the matter power spectrum 
is determined mostly by one parameter, k eq = a eq H eq oc Q m /z 2 . 
The matter power spectrum also scales with the primordial am¬ 
plitude ,4 S ; keeping A s fixed, but increasing k eq , means that the 
entire spectrum shifts sideways so that lenses of the same typ¬ 
ical potential depth v Fi ens become smaller. Theoretical ACDM 
models that keep { eq = k eq x* fixed will therefore have the same 
number (proportional to k eq x *) of lenses of each depth along 
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Fig. 7 Parameter constraints from CMB lensing alone in the 
ACDM model (samples, colour-coded by the value of the 
Hubble constant) using the minimal priors described in the text; 
grey bands give the corresponding 1 cr and 2 cr constraints us¬ 
ing the approximate fit of Eq. (13). Solid coloured contours 
show 68 % and 95 % constraints when additional information 
is included: BAO data from SDSS and 6DF (Anderson et al. 
2014; Ross et al. 2015; Beutler et al. 2011; blue), and the same 
but fixing the CMB acoustic-scale parameter (9\ic to a CMB 
power spectrum fit (red; (9\tc = 1.0408). Solid black contours 
show the constraint from the Planck CMB power spectra (Planck 
temperature-! low-7 polarization). 


the line of sight, and distant lenses of the same depth will also 
maintain the same angular correlation on the sky, so that the 
shape of the spectrum remains roughly constant. There is there¬ 
fore a shape and amplitude degeneracy where £ eq ~ constant, 
A s ~ constant, up to corrections from sub-dominant changes in 
the detailed lensing geometry, changes from late-time potential 
decay once dark energy becomes important, and nonlinear ef¬ 
fects. In terms of standard ACDM parameters around the best-fit 
model, Cq “c with the power-law dependence on O m only 

varying slowly with Q m ; the constraint ( eq oc Q^ 6 /t = constant 
defines the main dependence of Hq on Q m seen in Fig. 7. 

The argument above for the parameter dependence of the 
lensing power spectrum ignores the effect of baryon suppres¬ 
sion on the small-scale amplitude of the matter power spectrum 
(e.g. Eisenstein & Hu 1998). As discussed in Appendix E, this 
introduces an explicit dependence of the lensing power spectrum 
on Q. m h 2 . However, since the parameter dependence of f eq is 
close to (Q m h 2 ) 1/2 , we can still think of A s and ( eq as giving 
the dominant dependence of on parameters. 

In practice the shape of C'^ is not measured perfectly, and 
there is also some degeneracy between A s , which directly in¬ 
creases the amplitude, and £ eq , which increases the amplitude via 
an increase in the number of lenses along the line of sight and in¬ 
directly by changing the depth of potential wells of a given size. 
On small angular scales with L > £ cq , 'F| ens becomes a strong 
function of scale due to the fall-off in the matter power spectrum 
for k > k eq . For a fixed comoving scale, the amplitude increases 
if k eq becomes larger since this reduces the amount of decay of 
the gravitational potential during radiation domination (Pan et al. 
2014), and also lessens the impact of baryon suppression. This 


increases the dependence of C, on £ eq beyond the linear scaling 
from the number of lenses, and on small scales L 4 C^ oc A s C\ q " L , 
where ni > 0 determines the actual /,-dcpcndcnce (inherited 
from the scale dependence of the small-scale matter power spec¬ 
trum and the effect of baryon suppression). For Planck, the lens¬ 
ing spectrum is best measured at L ~ 200: since «/ increases 
with L, making A s larger and t). q n20 ° smaller to keep the power 
constant at L ~ 200 leads to an increase in power at lower L and 
a decrease in power at higher L (see Fig. 8). The lensing ampli¬ 
tude at the peak L « 50 (and hence the rms deflection angle), is 
therefore slightly anti-correlated with the amount of small-scale 
power (as shown by the lensing only constraint in Fig. 9; see 
discussion below). 

The actual three-parameter constraint is close to 

A s oc A s 7 2 q 3 oc constant + 3 %, (14) 

which, in terms of A s , £ cq , and h is independent of h, as expected, 
given the parameter dependence of discussed above. Recall 
that the fairly high power-law dependence with £ eq comes from 
the rapid fall-off in the matter power spectrum with scale (and 
from the dependence on Q. m h 2 of baryon suppression effects). 
Approximately the same scaling will also affect other measures 
of small-scale power, for example erg. From Eq. (E.20), for pa¬ 
rameter values of interest erg and A s are related approximately 
by 

^Q^ckA^ 5 . (15) 

The parameter direction crgfij]; 25 = constant is therefore close 
to the direction in which A s (and hence L 4 C^) is constant for 
fixed 4q- The combination A s £l q is also not that far from A s £ 2 q 
of Eq. (14), which remains nearly constant as ( cq and A s vary 
within their degeneracy. The crgQ,^ 4 degeneracy direction shown 
in Fig. 7 and the constraint of Eq. (13) are therefore significantly 
tighter than the corresponding constraint in the A s -Q m plane, be¬ 
cause erg absorbs some of the degeneracy between matter density 
and primordial power at a fixed scale. 

In terms of <xg, a well-determined three-parameter combina¬ 
tion is 1 

crgQf n 25 (Q m h 2 Y° 31 = 1.228 ± 0.028 (lensing only; 68 %), 

(16) 

at 2.3% precision. In terms of power, the 4.6% precision of 
the square of this parameter combination is still weaker than 
the template amplitude measurement of Eq. (11a). This is ex¬ 
pected, given the approximation of a power-law fit, and that in 
Eq. (16) we have also marginalized over variations in n s and 
Qb h 2 . Marginalizing out the Q m /r dependence in Eq. (16), the 
constraint on crgQ^ 25 of Eq. (13) is further degraded to 3.5 %. 

We can also quantify the amplitude of lensing by the rms 
deflection angle ( d 2 )^ 2 , which we define via 

2000 nr m n 

(d 2 ) = £ + 1)C l ■ (17 > 

L —2 

This is approximately an integral over all scales. It is propor¬ 
tional to the primordial fluctuation power A s , and scales roughly 
as a power law in £ eq . A simple argument based on the num¬ 
ber of lenses along the line of sight would suggest (d 2 ) oc A s 4q- 

4 This direction differs from that in Eq. (14), since the orthogonal 

direction in the A s -f eq plane is not that much worse constrained, so the 
Fisher information in the three-parameter space is not fully dominated 
by a single direction. 


9 







Planck Collaboration: Gravitational lensing by large-scale structures with Planck 



L 


■ 


0.152 
0.144 
0.136 
0.128 a 

n 

0.120 =£> 
0.112 
0.104 
0.096 



0.125 
0.124 
0.123 
0.122 
0.121 p 
0.120 ^ 
0.119 
0.118 
0.117 
0.116 


Fig. 8 Lensing potential power spectra in the base ACDM model drawn from the posterior of the lensing-only likelihood (left) and 
the Planck TT+lowP likelihood (right), coloured by the value of the dark matter density Q c /t 2 . For fixed CMB power spectra (i.e. the 
right-hand plot) the lensing power spectrum shape is nearly fixed, but the amplitude increases as the matter density increases across 
the acoustic-scale degeneracy (the accompanying reduction in A s is sub-dominant). With lensing reconstruction data, as Q. c lr varies 
the lensing spectrum changes shape, with the amplitude around the best-measured L ~ 150, remaining pinned by the accuracy of 
the measurement there. A combination of the lensing and Planck TT+lowP data measures the matter density significantly better, 
with the higher values in the right-hand plot being excluded. Note that the colour scales on the plots are different. 


This is almost correct, in particular the strong dependence of Cj 
on f eq integrates out when forming (d 2 ) (see Eq. E.4). However, 
there is a further dependence on O m and h due to the effects of 
baryon suppression and the late-time decay of the gravitational 
potential, which is important at low-L; we find approximately 

(d 2 )°cA s i: eq (nji) l/5 , ( 18 ) 

by taking derivatives of C^' computed with CAMB. The rms is 
measured at 2.4 % precision to be 

(d 2 ) 1 ^ 2 = (2.46 ± 0.06) arcmin (lensing only; 68 %). (19) 

So far we have only considered a ACDM model. However, 
for most generalizations that are close to ACDM, the same scal¬ 
ing also approximately holds, and the constraint direction of 
cr x £l||; 2S is only weakly model-dependent, although the centroid 
can shift and also the relation between Q m and h. For exam¬ 
ple, for a model with massless sterile neutrinos parameterized 
by (V e ff, the constraint of Eq. (13) is virtually identical, but pa¬ 
rameter samples shift towards systematically higher O m /z 2 than 
in ACDM to keep the equality scale roughly constant when 
N e ff > 3.046. As a further example, in models with three massive 
(active) neutrinos with degenerate masses, marginalizing over 
the mass of massive neutrinos shifts the constraint slightly down, 
giving 

cr 8 O^ 25 = 0.566 ± 0.025 (ACDM+Im v , lensing only; 68 %). 

( 20 ) 

Note that in the neutrino case non-linear corrections from 
halofit are unreliable away from Q,„ ~ 0.3, which can affect 
this result at around the 0.01 level. 

3.5.2. Joint parameter constraints 

The small-scale CMB power spectra are only directly sensi¬ 
tive to the primordial power A s after damping by reionization, 
through the combination A s e~ 2T . However lensing also smooths 
the power spectra at the several percent level, which allows the 


amplitude A s to be constrained via the amplitude of the smooth¬ 
ing effect (and a small transfer of power to small scales), even 
when r is not constrained precisely. 1 The lensing smoothing can 
be thought of as approximately a convolution of the deflection 
angle power spectrum with the CMB power spectrum, and hence 
is mostly sensitive to scales around the peak of the deflection 
angle power spectrum. We can quantify the amplitude of the 
smoothing approximately by the mean-squared deflection angle 
(d 2 ) from Eq. (17). 

The CMB power spectra determine the angular acoustic 
scale 0* * very accurately, which can be used to eliminate the h- 
dependence from (d 2 ), considered as a function <r x , h, and Q m . 
Empirically, we find that in the ACDM model the CMB power 
spectrum posterior has (d 2 ) 1 ^ 2 oc ct 8 D|J; 2 , close to the crgD^ 25 
combination measured more directly by lensing reconstruction; 5 
see Fig. 9. For a flat ACDM model the constraint on the latter 
parameter is 

cr 8 Q“ 5 = 0.621 ± 0.013 (Planck TT+lowP; 68 %). (21) 

The amplitude information here is coming from a combination 
of the lensing smoothing effect, and information from the CMB 
power spectrum amplitude combined with the constraint on r 
from low-7 polarization. The lensing-only constraint of Eq. (15) 
is slightly weaker, but somewhat lower than this (consistent at 
just over 1 cr); see Fig. 7. The combined constraint is tighter and 
restricted to slightly lower values compared to the Planck power 


^ The smoothing effect of lensing in the Planck TT power spectrum 
is detected at around 10 cr (Planck Collaboration XIII 2016). It is also 
detected at lower significance in the polarization power spectra. 

* The empirical parameter dependence of (d 2 ) 112 on for samples 
from the CMB power spectrum posterior differs from what one would 
infer from Eqs. (15) and (18). These equations assume fixed n s and 
kl b h 2 , which is appropriate for our lensing-only constraints with the as¬ 
sumed priors, but for the CMB samples the fractional dispersion in n s 
and £f b is comparable to that in cr 8 , h, and £f m and so cannot be ne¬ 
glected. 
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Fig. 9 Constraints on the rms lensing deflection angle (d 1 )^ 1 
from Planck TT+lowP CMB power spectra in the base ACDM 
model (samples, colour coded by the value of the Hubble con¬ 
stant), showing that the inferred lens deflection amplitude is 
strongly correlated with the parameter combination crgflj}; 25 . The 
solid black contour is the constraint from CMB lensing only, 
which prefers lower values of the lensing amplitude, giving the 
joint constraint shown by the red contours. Since the lensing 
reconstruction mainly constrains scales smaller than the peak 
of the deflection angle power spectrum, there is an additional 
degeneracy that makes (d 2 ) 1 ^ 2 slightly anti-correlated with the 
measurement of lensing power on smaller scales, which is ap¬ 
proximately proportional to (the square of) crgfi^ 25 

spectrum measurement: 

crgQ^ 25 = 0.609 ± 0.007 ( Planck TT+lowP+lensing; 68 %). 

( 22 ) 

Here, for joint constraints we no longer apply the additional pri¬ 
ors used in the lensing-only analysis, and allow for the full lin¬ 
earized CMB power spectrum dependence of the lensing renor¬ 
malization in the likelihood. The overlap region in the £2 m -crg 
plane also excludes some of the lower Hubble constant values 
allowed by the power spectrum data, giving a slight shift of the 
mean in the direction of the higher values preferred by some lo¬ 
cal measurements (Riess et al. 2011; Humphreys et al. 2013; 
Efstathiou 2014), as shown in Fig. 9. This could be an addi¬ 
tional indication that the lower Hubble constant values from the 
CMB spectra are at least partly a random statistical fluctuation, 
in which case joint constraints should be usefully closer to the 
truth. The joint constraint gives 

H 0 = (67.8 ± 0.9) km s _1 Mpc' 1 

(Planck TT+lowP+lensing; 68 %), (23) 

about 0.5 cr higher than Planck TT+lowP alone, and in good 
agreement with the constraint from Planck TT+lowP+BAO 
data. 

The measurement of the fluctuation amplitude from lensing, 
combined with the amplitude of the CMB power spectrum on 
intermediate and small scales, can also be used to constrain the 
reionization optical depth independently of low-7 polarization. 
The CMB power spectra themselves provide a weak constraint 
via the lensing smoothing, which is significantly improved by 
using the additional information in the lensing reconstruction. 


Figure 10 shows the constraint when Planck CMB tempera¬ 
ture power spectrum constraints, without low-C polarization, are 
combined with the lensing likelihood. We find 

t = 0.070 ± 0.024 (Planck TT+lensing; 68 %), (24) 

corresponding to a detection of a non-zero optical depth at more 
than 2cr, and a corresponding reionization redshift z re = 9.0 1 22 j 
(68%). These results are consistent with the baseline Planck 
power spectrum constraints including low-7 polarization, which 
give t = 0.078 ± 0.019 (Planck TT+lowP; 68%). Both results 
indicate a consistent downward shift in mean compared to the 
central value from WMAP (r = 0.089 ± 0.014; Hinshaw et al. 
2013), in agreement with lower predicted values based on recent 
results for the integrated luminosity function (e.g. Finkelstein 
et al. 2015). There is a corresponding downward shift in cr 8 , with 
joint constraint 

crg = 0.815 ± 0.009 (Planck TT+lowP+lensing; 68 %). (25) 
For further discussion see Planck Collaboration XIII (2016). 



Fig. 10 Constraints on the reionization (mid-point) redshift z re 
and erg from a combination of Planck CMB power spectra and 
the lensing reconstruction, excluding low-7 polarization infor¬ 
mation. The grey band denotes the approximate z re < 6.5 re¬ 
gion excluded by observations of the spectra of high-redshift 
quasars (Fan et al. 2006). The lensing data significantly shrink 
the allowed region of parameter space, preferring lower values 
of the reionization redshift, in good agreement with the measure¬ 
ment from Planck low-7 polarization. 


3.5.3. Constraints on extensions to the ACDM model 

The CMB power spectra only weakly constrain extended mod¬ 
els that provide multiple ways to combine late-time parame¬ 
ters to give the same observed angular diameter distance to 
recombination. The main Planck CMB power spectrum con¬ 
straints on these parameters are driven largely by the impact of 
lensing smoothing on the acoustic peaks, which weakly breaks 
the geometric degeneracy. As discussed in detail in Planck 
Collaboration XIII (2016), the power spectrum results appear to 
prefer larger lensing smoothing than in ACDM (in this release. 
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Fig. 11 Constraints on the curvature and rms lensing deflec¬ 
tion angle from Planck TT+lowP CMB power spectra (samples, 
colour coded by the value of the Hubble constant), showing that 
the closed-universe models favoured at the 2 cr level have high 
lensing power that is inconsistent with the lensing reconstruc¬ 
tion result. The combined lensing result (black contours) is much 
tighter, and provides additional information on top of the con¬ 
straint when adding BAO data (red contours); both joint con¬ 
straints, and the combined filled contours, are consistent with a 
flat universe. 


slightly above the 2 cr level). It is therefore particularly useful to 
study joint constraints with the lensing likelihood, which con¬ 
strains the lensing amplitude more directly and more tightly. 
Since the lensing-reconstruction power spectrum has an ampli¬ 
tude that peaks slightly lower than the ACDM predictions from 
the CMB power spectra, the joint constraint eliminates a large 
region of parameter space with positive curvature allowed by 
the power spectrum results. For example. Fig. 11 shows the con¬ 
straints on the curvature parameter O-k in non-flat models; the 
high lensing amplitudes, as show in the figure by the rms de¬ 
flection angle, are ruled out by the lensing reconstruction, and 
the joint constraint is consistent with a flat universe {D-k = 0). 
The improvement in error bars is dramatic, with the CMB power 
spectra giving 

Q k = -0.052!°;“ (Planck TT+lowP; 68 %), (26) 

and the joint constraint 

Q. k = -0.005!°;°°° C Planck TT+lowP-t-lensing; 68 %). (27) 

This measurement of curvature from the CMB alone has sub¬ 
percent precision. External data can also break the geometric 
degeneracy; the joint BAO contours are shown in Fig. 11 and 
are consistent. Since the BAO results are independent of fluctu¬ 
ation amplitude, the lensing-reconstruction joint constraint adds 
significant additional information, shrinking the allowed param¬ 
eter space compared to Planck +BAO alone. 

In addition to curvature, CMB results for dark energy mod¬ 
els also suffer a geometric degeneracy, and lensing results can 
similarly provide additional information; for applications to a 
variety of dark energy models and modified gravity see Planck 
Collaboration XIV (2016). Lensing can also improve the limits 
on many other parameters, though in a less dramatic fashion if 


the CMB power spectra are themselves able to place a good con¬ 
straint. Full results for many combinations are available online 
in the grid tables and parameter chains. 1 

3.5.4. Sensitivity to L-range 

The results in the previous subsections were obtained using our 
conservative baseline likelihood, with eight bins in the multi¬ 
pole range 40 < L < 400. We now briefly discuss how joint 
constraints from the Planck power spectra and lensing change 
for different cuts of the lensing data. 

Using the first nine bins of the aggressive binning de¬ 
scribed in Table 1, covering 8 < L < 425, results from 
Planck TT+lowP+lensing are consistent at the 0.1 cr level with 
the conservative likelihood, in both ACDM and ACDM+ £ m v 
cosmologies. Extending the range to higher L, the Planck 
TT+lensing constraint on the optical depth is stable, however 
the measured amplitude <x 8 fl[! 25 shifts downwards by about 1 cr 
with a negligible change in the uncertainty. About half of this 
comes from including the outlying bin at 638 < L < 762. In the 
ACDM+2 m v model, this preference for lower amplitudes on 
smaller scales pulls joint neutrino mass constraints up by about 
1 cr. The baseline result is 

^ m v < 0.145 eV (Planck TT+lowP+lensing+BAO; 68 %), 

^ (28) 
but using the full L-range the posterior peaks slightly (though 
not significantly) away from zero: 

^ m v = 0.16!o'°j eV (Planck TT+lowP+aggressive lensing 

+BAO; 68 %). 

(29) 

As discussed in Sect. 4, the multipole range 300 < L < 900 may 
be unreliable, specifically due the failure of the TTTT curl test. 
A priority for future analysis should be to develop a trustworthy 
likelihood over the full multipole range. 


4. Consistency and null tests 

With multiple frequency bands and the combination of both tem¬ 
perature and polarization data, the Planck full-mission data set 
provides many opportunities to test internally the lens recon¬ 
struction. In this section we present consistency and null tests 
for the MV lensing potential presented in Sect. 3. 

For consistency tests, we perform lens reconstruction using 
different data/analysis choices than those used for our fiducial 
MV reconstruction, obtaining an alternative measurement of the 
lensing potential power spectrum C^j test , which we then com¬ 
pare to the baseline MV reconstruction. To compare the test and 
MV reconstructions more quantitatively, we use a^ 2 statistic for 
the difference between the bandpowers, which is calculated as 


Vtcst 


S' 




L I test 


C' 


L I MV 


s'f (d*f\ 


A(p(p\ \ 

^L' ImV/ ’ 


(30) 

where £ tesl is a diagonal covariance matrix for the difference 
obtained from simulations of the two lens reconstructions. We 
compare this x 1 to a set of values determined from simulations, 
to estimate a PTE. We do this for several possible choices of 
binning: a single bin for either the conservative (40 < L < 400) 
likelihood range, or the low-L and high-L ranges (8 < L < 40 


t http://www.cosmos.esa.int/web/planck/pla 
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and 400 < L < 2048, respectively) that are additionally included 
in the aggressive likelihood; and the individual bins that span the 
conservative and high-L multipole ranges. The single-bin tests 
check for broad-band differences between the two reconstruc¬ 
tions, while the multiple-bin tests provide a more sensitive check 
on shape differences or inaccuracies in the error bars for the two 
reconstructions. 

We group the consistency tests as follows, based on the is¬ 
sues which they are designed to probe. 

• ‘Individual estimator crosses’ that test for consistency 
among the five quadratic estimators used to measure (f>. 

• ‘Foregrounds tests’ that use single-frequency maps or differ¬ 
ent component-separation methods to test for possible fore¬ 
ground contamination. 

• ‘Noise tests’ that check that our power spectrum estimates 
are not biased by instrumental noise, using cross-estimators 
designed to avoid noise biases. 

• ‘Analysis choices’ that test our sensitivity to choices such as 
the unmasked sky fraction and multipole range. 

• ‘Bias-hardened estimators’ in which we use lensing estima¬ 
tors that should be less sensitive to certain systematic and 
instrumental effects. 

In Fig. 12 we plot several of the most important consistency 
tests, as well as their difference from the MV reconstruction. 
The^ 2 values associated with the full set of consistency tests are 
presented in Table 2. The tests themselves are described in more 
detail in Sects. 4.1-4.5. 

In addition to consistency tests, we also perform several null 
tests in which we apply our lens reconstruction and power spec¬ 
trum estimation procedure to maps that have been differenced 
to remove sky signal, or use curl-mode lensing estimators that 
should have zero contribution from the gradient lensing poten¬ 
tial. We check that the overall amplitude (for a template given by 
the fiducial theory spectrum C* /, Md ) is consistent with zero, and 
also evaluate a^f 2 computed as 

xU =S?(<21 (31) 

where E nu u is a diagonal covariance matrix measured from simu¬ 
lations of the particular null test. A PTE for each null test is com¬ 
puted by comparison to simulations in the same way as for the 
consistency tests. The null power spectra are plotted in Fig. 13, 
and summary statistics are given in Table 3. The null tests are 
discussed in more detail in Sect. 4.6. 

4.1. Individual estimator crosses 

With the foreground-cleaned temperature and polarization maps 
from SMICA we can form five distinct lensing estimators, and 15 
distinct auto- and cross-spectra. We plot all of these spectra indi¬ 
vidually in Fig. 14. The lensing amplitudes, as well as PTEs for 
the differences from the MV reconstruction are given in Table 2. 
There is a wide variation in the sensitivity of the different es¬ 
timators, although we see high (greater than 5 cr significance) 
detections of lensing in nine of the spectra. ’ 

: In the 2015 Planck data release, there is some evidence in the polar¬ 
ization power spectra for uncorrected systematic effects associated with 
temperature-to-polarization leakage (Planck Collaboration XIII 2016; 
Planck Collaboration XI 2016). Although we find no evidence for in¬ 
consistencies in the polarization-based lens reconstructions, the statis¬ 
tical errors are large and may obscure small systematic effects that are 
subdominant for the current analysis. Leakage from the CMB temper- 


As a further test, we combine several of the spectra with the 
highest S/N that only use estimator crosses (TTxEE, TTxTE, 
TT x EB, and EE xTE) into a single ‘crosses’ estimator. This 
combination of cross-spectra has an order-of-magnitude smaller 
disconnected noise bias correction than the auto-spectrum of 
the fiducial MV estimator, with AA™’ 0 Q 4 g ises = 0.4, compared to 
A/4^° p ^ = 4.3 for the baseline MV lensing reconstruction. The 
results of this test are presented in Fig. 12 and Table 2. With 
the crosses power spectrum estimates we are able to test that 
the lensing amplitude is not biased by errors in the disconnected 

noise bias correction AC^Ino at around the 2% level. 

We also form a combined power spectrum estimate using 
only polarization data (the EE and EB estimators). This ‘Pol. 
only’ estimator gives an amplitude estimate for the conservative 
multipole range of A = 0.76 ± 0.15, a 5 cr detection of lensing at 
a level consistent with the fiducial MV reconstruction. 

4.2. Foreground tests 

The SMICA CMB map that is used for our baseline MV recon¬ 
struction combines multiple frequency bands as a function of 
angular scale to reject foregrounds and to produce a single high- 
fidelity estimate of the CMB. We can also compare to a simpler 
approach based on individual frequency maps, which we clean 
by subtracting a dust template. Our baseline results in Planck 
Collaboration XVII (2014) were based on this approach using 
143 GHz and 217 GHz data, with the 857 GHz map projected out 
during the filtering process as a dust template. In this work we 
repeat this analysis. We use 857 GHz as a dust template in tem¬ 
perature, and no cleaning in polarization. The results of this test 
are presented in Fig. 12 and Table 2 with the labels ‘143 GHz’ 
and ‘217 GHz’, where it can be seen that we obtain lensing band- 
powers that are consistent with the MV reconstruction for the 
conservative multipole range. In the simulations used to char¬ 
acterize these reconstructions, we add coloured Gaussian noise 
to the 143 GHz and 217 GHz temperature map simulations, de¬ 
signed to mimic the power of unsubtracted extragalactic fore¬ 
grounds. For this we use the same source model used in Planck 
Collaboration XVII (2014). We also add Gaussian power to the 
polarization maps, designed to mimic the power due to Galactic 
dust contamination that can be seen on large angular scales. The 
resulting temperature and polarization simulations have power 
spectra that agree with those of the data at the percent level in 
power. 

We can also test the stability of the lensing reconstruction 
to the choice of the SMICA component-separation method over 
other alternatives such as NILC, SEVEM, and Commander. The 
results of these comparisons are also shown in Table 2; all are 

ature anisotropies could impact the lensing analysis in several ways. 
First, the disconnected (Gaussian) part of the CMB 4-point function 
is modified, and this would affect both our mean-field and NO correc¬ 
tions (see Appendices A.2 and A.3). These corrections are made us¬ 
ing a mix of data and simulations, such that even if the simulations 
do not fully capture the leakage (note they do include leakage due 
to main-beam mismatch, for example), the additive bias to the lens¬ 
ing power spectrum is only second-order in any unmodelled effects. 
The second way that leakage can affect the lensing power spectrum is 
through changes to the lensing response functions in Eq. (A. 13), which 
propagate through to normalization errors. This should be partially cap¬ 
tured by a Monte-Carlo correction that we make to the lensing power 
spectrum. We have not attempted to quantify more carefully any resid¬ 
ual effects of temperature-to-polarization leakage, since, ultimately, the 
polarization data carries relatively little weight in our final MV recon¬ 
struction. 
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Fig. 12 Summary of lensing power spectrum consistency tests. Upper panel: comparison of the lensing power spectra for the 
baseline MV reconstruction and several alternative estimators with different data/analysis choices discussed in the text. Lower 
panels: detailed plots of the difference between each alternative reconstruction and the MV power spectrum, in units of the MV 
error bars (cr). The ±lcr band is indicated in grey. Error bars on individual data points give the standard deviation of the expected 
difference from the MV reconstruction for each case, estimated using Monte Carlo simulations. 
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Fig. 13 Power spectrum null tests discussed in Sect. 4.6. Curl null tests are denoted by if/. The individual null test power spectra 
have been scaled by the factors indicated in each panel for visual purposes. The fiducial (gradient-mode) lensing power spectrum is 
shown as the dashed line in each panel. 
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Fig. 14 Grid of lensing potential power spectra for the 15 possible auto- and cross-spectra of the five quadratic lensing estimators 
obtained from the SMICA foreground-cleaned maps. The measurements are arranged in order of their total S/N for a measurement 
of the lensing amplitude. The MV power spectrum estimate, which combines the individual cross-spectra and is the basis for our 
fiducial lensing likelihood, is plotted in the top-left panel. 


very consistent with the SMICA-based reconstruction over the 
conservative multipole range, although the PTEs for the over¬ 
all amplitude at high-L (400 < L < 2048) are low for NILC 
and Commander. For the SEVEM and Commander analyses, po¬ 
larization maps were not available at the high resolution used 


here ((V S ide = 2048), and so we have limited the comparisons to 
temperature only. 

As a test of Galactic contamination, we also repeat our base¬ 
line analysis (performed on a fraction /) ky = 0.7 of the sky) 
using a more aggressive (/A y = 0.8) and more conservative 
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(/sky = 0.6) Galactic mask. These masks are constructed by 
thresholding maps from the Planck high-frequency channels to 
achieve the desired sky fraction. In Table 2 we see that for both 
alternative sky fractions we estimate lensing bandpowers that are 
consistent with the baseline. 

As a final and very powerful test of high-frequency fore¬ 
ground contamination, we also perform a lens reconstruction in 
which we use the 353 GHz map as one of the four legs in the 
trispectrum used to estimate C' 1 " 1 ’. Based on the approximate fre¬ 
quency scaling of dust (Gispert et al. 2000) we expect that any 
dust foreground bias 7? ( ] llst ’ to the lensing power spectrum esti¬ 
mate that we see in this test will scale roughly as 


b 


353 

dust 


lh 217 


* 20 b 


143 

dust’ 


(32) 


and so this should provide a sensitive test for such contami¬ 
nation in the SMICA map. The 353 GHz map has significantly 
more power than CMB+noise-only simulations, due primarily 
to Galactic dust, and temperature fluctuations from the CIB. To 
obtain reasonable error bars from the 353 GHz simulations used 
to characterize this reconstruction, as with the other reconstruc¬ 
tions we add an isotropic Gaussian sky signal designed to match 
the simulation power spectrum to a smoothed power spectrum of 
the data. When filtering the 353 GHz map, we use the same pa¬ 
rameters as used for the SMICA map (except for the beam trans¬ 
fer function), to have the same weight given to different angular 
scales in the quadratic lensing estimators. As with the 143 GHz 
and 217 GHz analyses, we project out the 857 GHz map as a 
dust template from the 353 GHz map in temperature, but make 
no attempt to subtract foreground contamination in polarization. 
As can be seen in Fig. 12 and Table 2, we obtain good agree¬ 
ment between the ‘one leg 353 GHz’ reconstruction and the MV 
lensing power spectrum. Lensing power is detected at 20 cr even 
using 353 GHz as one of the four estimator legs. The stability of 
the lens reconstruction in the presence of the foreground power 
at 353 GHz is very striking. As a further test for the stability 
of the polarization-based lens reconstructions, we can also use 
the ‘Pol. only’ estimator from Sect. 4.1 with one leg given by 
the 353 GHz data. This reconstruction is very noisy, however we 
do obtain an overall amplitude as well as bandpowers that are 
consistent with the MV ‘Pol. only’ reconstruction, as shown in 
Table 2. 


4.3. Noise tests 

As mentioned in Sect. 2, the FFP8 simulations on which our 
results are based underestimate the noise power by several per¬ 
cent at high (. We have dealt with this problem in our simu¬ 
lations by including isotropic Gaussian power to make up for 
the difference. This neglects the fact that the missing noise 
component most likely has an anisotropic distribution (which 
could lead to errors in the mean-field subtraction), and could 
be non-Gaussian as well. To test that our treatment of the miss¬ 
ing noise component is adequate, we have performed several 
tests. Most stringently, we form quadratic lensing estimates 
using cross-correlation between pairs of maps with indepen¬ 
dent noise realizations. We constmct both a frequency cross 
143 GHz x 217 GHz, as well as a cross between the first and sec¬ 
ond halves of the mission data. These estimators have no noise 
mean field, and the lensing power spectmm estimated from them 
is insensitive to non-Gaussianity of the instrument noise. The re¬ 
sults of this analysis are given in Table 2, where we see results 
that are consistent with the MV power spectrum estimate. 

As an additional test, we have recreated the analysis of 
Planck Collaboration XVII (2014), where the missing noise 


power was accounted for by scaling the noise component of the 
FFP8 simulations by a small factor. This approach reasonably 
supposes that the anisotropy of the missing noise component 
is the same as that of the modeled component. This test is la¬ 
belled ‘Noise scaling’ in Table 2. Again, we see a power spec¬ 
trum measurement with an amplitude that is very similar to the 
MV estimate (although we note that because the only difference 
in this test is the constmction of the simulations that are used, 
the expected scatter between the ‘Noise scaling’ and baseline 
approaches is difficult to determine). 

4.4. Analysis choices 

We test the stability of our reconstruction to several analysis 
choices, re-running the baseline analysis implementing the fol¬ 
lowing changes one at a time. 

• Use a high-pass filter in temperature and polarization with 
Aniri = 1000 rather than £ mm = 100, as adopted in the baseline 
analysis. 

• Use a low-pass filter with f’ max = 1500 rather than 
/’max = 2048, as used in the baseline analysis. 

• Apodize the lensing convergence estimate before power 
spectrum estimation, following the analysis in Planck 
Collaboration XVII (2014). This change reduces the size of 
the Monte Carlo correction, but does not produce a signifi¬ 
cant change in the lensing bandpower estimates. 

The results of these tests are given in Table 2. 

4.5. Bias-hardened estimators 

Throughout this paper, we have used quadratic estimators opti¬ 
mized to detect the anisotropic covariance between CMB modes 
that is induced by (fixed) lensing. These estimators are biased 
by non-lensing sources of such mode coupling, including ef¬ 
fects such as galaxy and point-source masking, inhomogeneous 
instrumental noise, and beam asymmetries. As discussed in 
Sect. 2, we estimate these biases using Monte Carlo simula¬ 
tions and subtract them as a ‘mean-field’ contribution from the 
estimated lensing potential (as in Eq. A. 17). A complementary 
method, which is less reliant on the quality of the Monte Carlo 
simulations used to perform this debiasing, is the bias-hardening 
procedure advocated by Namikawa et al. (2013). In this ap¬ 
proach, one constmcts a new quadratic lensing estimator that is 
orthogonalized to worrisome sources of bias. The constmction 
of these estimators is reviewed in Appendix B. In the analysis 
here we focus on temperature-only bias-hardened estimators to 
investigate the effects of bias-hardening, because mean-field cor¬ 
rections are much smaller in polarization than in temperature. 
We denote these estimators as rf> TT - X, where X denotes the 
bias(es) that are orthogonalized against. Here, we consider the 
biases due to masking (M), inhomogeneous noise ( N ), and the 
shot noise due to unresolved point sources ( S ). As discussed in 
Osborne et al. (2014), the source-hardened estimator <f> TT -S also 
reduces contamination from effects such as the cross-correlation 
between the unresolved sources and the lensing potential. 

Mean-field corrections are largest at low-L (the power spec¬ 
trum of the mean-field correction is larger than the fiducial lens¬ 
ing potential power spectmm for L ^ 100), and so we be¬ 
gin by looking closely at these scales. In Fig. 4.5 we plot the 
power spectra of the mean-fields for each of the bias-hardened 
estimators. In the lower panel we also plot the difference be¬ 
tween each estimator and the non-bias-hardened TT estimator. 
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Fig. 15 Upper, power spectrum of the low-L mean-field (mea¬ 
sured from Monte Carlo simulations) for the TT lensing esti¬ 
mator as well as several bias-hardened estimators. The mean- 
field power spectrum of the TE lensing estimator is plotted 
for comparison, and the fiducial lensing potential power spec¬ 
trum is shown in dashed black. To avoid visual clutter we have 
not shown the source-hardened estimator TT - S ; its spectrum 
is midway between the noise-hardened (TT - N) and mask- 
hardened (TT - M ) cases. Lower, power spectra of the differ¬ 
ences between the TT lensing estimate and the bias-hardened 
reconstructions or the TE reconstruction (with the same colour 
scheme as in the top panel). The expected power spectra of the 
differences, estimated from simulations, are plotted as black dot¬ 
ted lines. 


For L > 8, we see agreement between all of the bias-hardened 
estimators and the non-bias-hardened reconstruction, regardless 
of the much smaller mean-field corrections for the former. We 
have also used the bias-hardened estimators to measure the lens¬ 
ing potential power spectrum on smaller scales. In Fig. 12 and 
Table 2 we report lensing bandpower estimates using the bias- 
hardened estimators, which we find to be consistent with the 
baseline MV reconstruction for L < 400. 

4.6. Null tests 

In addition to the consistency tests described above, we have 
also performed several direct null tests, either using estimators 
or maps for which the lensing signal is expected to be zero. 
These are plotted in Fig. 13, and significance statistics for each 
are given in Table 3. We discuss them in more detail below. 


mately 85°. Each scan ring is observed for between 39 and 65 
minutes, and then the satellite is repointed. This leads to a nat¬ 
ural ‘half-ring’ null test in which we run our standard lensing 
estimator on the half-difference of maps constructed from the 
first and second half of each scanning ring, which should consist 
entirely of noise. The PTE for the spectrum of the lensing recon¬ 
struction estimated from these half-ring difference maps is rather 
low for the conservative multipole range (see Table 3). Since the 
maps should be only noise by construction, to obtain an accept¬ 
able x 2 for the power spectra of reconstructions based on such 
maps requires very good modeling of the instrument noise. As 
can be seen from the factor of 200 scaling applied in Fig. 13, 
the impact on the reconstructed of minor errors in the noise 
modeling that likely drive the low PTEs for the null test should 
be negligible. 


Half-mission difference: Noise that is correlated between rings 
(for example, due to glitches and destriping) can be nulled in 
the half-ring difference. Another differencing approach, which 
avoids this problem, is to divide the entire mission into two 
halves that are then differenced. This half-mission difference 
requires a slightly larger analysis mask, to account for the fact 
that the sky coverage is imperfect in both the first and second 
halves of the mission and some pixels are missed. The PTEs for 
the half-mission test are reasonable. 

Curl-mode estimators: In general, deflections due to lensing 
can be decomposed into gradient ((/>) and curl (if/) modes. The 
curl modes are expected to be consistent with zero for lens¬ 
ing by scalar perturbations at Planck reconstruction noise levels. 
We construct curl-mode power spectrum estimators, analogous 
to those used for the gradient lensing modes, and measure their 
power spectra. Full-sky weight functions for these estimators are 
presented in Namikawa et al. (2012). Note that although they do 
not couple directly to the lensing potential, the curl power spec¬ 
tra do receive ‘Nl’-type contributions from the lensing trispec¬ 
trum, which we estimate analytically and subtract as with the 
Cf power spectrum estimates. We do not make ‘PS’ or ‘MC’ 
corrections for the curl-mode power spectra. In Table 3 we re¬ 
port estimated amplitudes for a signal with the shape of the fidu¬ 
cial lensing power spectrum C^' fid in the curl estimates. These 
amplitudes are all consistent with zero at better than 2 cr, with 
the exception of the TT reconstruction over the high-L range 
(400 < L < 2048). Indeed, there is mild evidence for a broad fea¬ 
ture in the TT curl-mode reconstruction centred around L = 500 
(see Fig. 13), which drives the low PTE of the TT curl-mode 
power spectrum over the high-L multipole range. Considering 
TT alone, the average curl power is low by about 2.9 cr in the 
high-L range, which is suggestive but not clearly unacceptable 
given the number of tests performed. Combined with the polar¬ 
ization reconstructions, which on their own appear consistent, 
the MV estimate remains low but only at 2 cr over the high-L 
range. Since the origin of this feature in the curl-mode spectrum 
is currently not understood, we have chosen to adopt the con¬ 
servative multipole range for the baseline Planck 2015 lensing 
likelihood. This multipole range has little overlap with that cov¬ 
ered by the curl-mode feature and was determined in advance of 
the curl result as it retains the majority of the S/N in the recon¬ 
struction. 


Half-ring difference: The Planck observation strategy consists 
of scanning the sky in rings with an opening angle of approxi- 
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Table 2 Amplitude fits for a suite of alternative reconstructions of the lensing power spectrum and comparison to the baseline MV lensing power spectrum estimate. The 
measured overall lensing amplitude A is given for the conservative ( 40 < L< 400) multipole range, as well as the lower signal-to-noise bins on either side. The difference between 
the test amplitudes and the corresponding MV result (with amplitudes of = 1.023 ± 0.085, A^^ 0() = 0.987 ± 0.025, and •A^JJ 204g = 0.788 ± 0.097 for the three multipole 
ranges) is given as Amv- The quoted uncertainty on the difference is obtained using Monte Carlo simulations and therefore accounts for any correlations between the MV and test 
estimates. The PTE values are estimated by comparison of the x 1 values to the distribution from simulations. PTEs are given for the square of the overall amplitude difference, 
as well as for the x 2 of the difference between the bandpowers of the binned spectra. 


Consistency test summary values 



A ± cr A 

Low-L 

8 < L < 40 

A M v ± °"a 

PTE Amv =o 

A ± 

: cr A 

Conservative 

40 < L < 400 

Amv ± ct a 

TTE Amv=0 

PTE bins 

A±ct a 

High-L 

400 < L < 2048 

Amv ± cr a 

PTE Amv=0 

PTEbi ns 

Individual Estimator Crosses 

TT x TT 

1.044 ± 

0.111 

-0.035 ± 0.075 

0.65 

0.990 

± 0.032 

-0.003 ± 0.023 

0.90 

0.47 

0.724 

± 0.096 

0.058 ± 0.042 

0.17 

0.81 

TT x TE 

1.000 ± 

0.108 

0.009 ± 0.070 

0.91 

0.969 

± 0.039 

0.018 ± 0.028 

0.54 

0.12 

0.628 

± 0.273 

0.155 ±0.243 

0.52 

0.52 

TTxEE 

0.872 ± 

0.149 

0.137 ±0.126 

0.29 

1.012 

± 0.055 

-0.025 ±0.051 

0.65 

0.89 

1.616 

± 0.448 

-0.833 ± 0.446 

0.06 

0.18 

TTxEB 

1.105 ± 

0.308 

-0.096 ± 0.295 

0.76 

0.988 

± 0.085 

-0.001 ± 0.081 

0.99 

0.83 

1.118 

± 0.270 

-0.335 ± 0.269 

0.22 

0.48 

TExTE 

0.890 ± 

0.194 

0.118 ±0.172 

0.47 

1.057 

± 0.079 

-0.070 ± 0.075 

0.35 

0.94 

3.607 

± 1.472 

-2.824 ± 1.461 

0.05 

0.41 

EExTE 

1.067 ± 

0.209 

-0.058 ±0.194 

0.77 

0.995 

± 0.085 

-0.008 ± 0.080 

0.93 

0.77 

-0.441 

± 1.961 

1.224 ± 1.964 

0.54 

0.18 

TExEB 

0.604 ± 

0.424 

0.405 ±0.415 

0.31 

1.016 

±0.157 

-0.029 ±0.156 

0.88 

0.86 

0.786 

± 1.235 

-0.003 ± 1.223 

1.00 

0.57 

TT x TB 

0.496 ± 

0.631 

0.512 ±0.630 

0.42 

0.890 

±0.180 

0.097 ± 0.178 

0.59 

0.76 

1.413 

± 0.575 

-0.630 ±0.571 

0.26 

0.06 

EExEE 

1.132 ± 

0.407 

-0.124 ±0.397 

0.72 

0.748 

±0.189 

0.239 ±0.186 

0.20 

0.50 

4.503 

± 5.723 

-3.720 ± 5.717 

0.49 

0.56 

EExEB 

1.530 ± 

0.664 

-0.522 ± 0.659 

0.44 

0.643 

± 0.247 

0.344 ± 0.246 

0.15 

0.20 

4.131 

± 2.432 

-3.348 ± 2.429 

0.19 

0.29 

TExTB 

2.426 ± 

0.837 

-1.417 ±0.833 

0.09 

0.722 

± 0.297 

0.266 ± 0.294 

0.37 

0.75 

5.242 

± 2.431 

-4.459 ± 2.425 

0.07 

0.20 

EExTB 

-0.981 ± 

1.212 

1.990 ± 1.209 

0.10 

1.288 

± 0.463 

-0.301 ± 0.463 

0.52 

0.01 

3.154 

± 4.492 

-2.371 ±4.482 

0.61 

0.81 

EBxEB 

0.588 ± 

1.819 

0.421 ± 1.817 

0.82 

1.401 

± 0.522 

-0.414 ± 0.521 

0.42 

0.39 

-0.422 

± 2.274 

1.205 ±2.263 

0.60 

0.80 

TBxEB 

0.887 ± 

2.954 

0.122 ±2.947 

0.95 

0.618 

± 0.798 

0.369 ± 0.798 

0.63 

0.34 

-2.591 

± 2.955 

3.374 ± 2.945 

0.26 

0.43 

TBxTB 

13.790 ± 

8.353 

-12.781 ±8.351 

0.12 

-2.809 

±2.189 

3.796 ±2.189 

0.09 

0.08 

-0.532 

± 8.268 

1.315 ±8.260 

0.86 

0.20 

Crosses 

0.958 ± 

0.097 

0.051 ±0.053 

0.33 

0.989 

± 0.033 

-0.002 ± 0.022 

0.91 

0.65 

1.028 

±0.193 

-0.245 ±0.174 

0.15 

0.84 

Pol. only 1 

1.252 ± 

0.350 

-0.243 ± 0.339 

0.46 

0.761 

±0.145 

0.226 ± 0.142 

0.11 

0.39 

1.787 

± 1.598 

-1.004 ± 1.586 

0.53 

0.84 

Foreground tests 

143 GHz 

0.958 ± 

0.096 

0.029 ± 0.069 

0.64 

1.024 

± 0.034 

-0.034 ± 0.024 

0.16 

0.22 

0.583 

±0.122 

0.198 ±0.093 

0.04 

0.02 

217 GHz 

0.941 ± 

0.092 

0.046 ± 0.077 

0.57 

0.963 

± 0.034 

0.027 ± 0.024 

0.30 

0.71 

0.929 

±0.109 

-0.148 ±0.088 

0.08 

0.80 

NILC 

0.977 ± 

0.073 

0.010 ±0.027 

0.75 

0.988 

± 0.026 

0.002 ±0.011 

0.83 

0.94 

0.877 

± 0.082 

-0.096 ±0.031 

0.01 

0.09 

sevem rr 11 

1.015 ± 

0.108 

0.016 ± 0.053 

0.78 

1.000 

± 0.032 

-0.011 ± 0.017 

0.53 

0.63 

0.758 

± 0.094 

-0.039 ± 0.047 

0.40 

0.73 

Commander 7T 11 

0.899 ± 

0.124 

0.132 ±0.092 

0.15 

1.018 

± 0.038 

-0.029 ± 0.028 

0.30 

0.42 

0.950 

±0.113 

-0.230 ± 0.099 

0.02 

0.21 

/sky = 0.8 

1.009 ± 

0.064 

-0.022 ± 0.028 

0.36 

0.987 

± 0.025 

0.003 ± 0.009 

0.73 

0.14 

0.761 

± 0.086 

0.019 ±0.031 

0.57 

0.22 

/sky = 0.6 

1.000 ± 

0.078 

-0.014 ± 0.034 

0.71 

0.991 

± 0.026 

-0.000 ± 0.010 

0.97 

0.90 

0.777 

±0.091 

0.004 ± 0.037 

0.92 

0.68 

(SMICA) a 353 GHz 

1.000 ± 

0.133 

-0.013 ±0.114 

0.93 

0.970 

± 0.048 

0.021 ± 0.041 

0.64 

0.76 

0.946 

±0.134 

-0.165 ±0.128 

0.21 

0.05 

(SMICA) 3 353 GHz pol. only 1 

0.277 ± 

0.996 

1.016 ± 1.076 

0.33 

0.264 

±0.451 

0.531 ±0.495 

0.27 

0.80 

7.045 

± 5.895 

-5.056 ± 6.440 

0.40 

0.31 

Noise tests 

143 x 217 GHz 

0.949 ± 

0.077 

0.038 ± 0.050 

0.44 

0.999 

± 0.027 

-0.008 ± 0.015 

0.61 

0.97 

0.755 

± 0.095 

0.025 ± 0.056 

0.66 

0.94 

HM1 x HM2 

1.028 ± 

0.081 

-0.041 ± 0.050 

0.39 

0.981 

± 0.029 

0.009 ± 0.016 

0.63 

0.67 

0.827 

± 0.090 

-0.047 ± 0.058 

0.40 

0.88 

Noise scaling 1 " 11 

0.989 ± 

0.070 

-0.002 ± 0.030 

0.93 

0.994 

± 0.026 

-0.003 ±0.011 

0.73 

1.00 

0.779 

± 0.084 

0.001 ± 0.037 

0.94 

1.00 

Analysis choices 

Cm = 1000 

1.136 ± 

0.104 

-0.149 ±0.081 

0.06 

0.990 

± 0.036 

0.000 ± 0.029 

1.00 

0.08 

0.718 

±0.153 

0.063 ±0.128 

0.57 

0.44 

Cax = 1500 

1.013 ± 

0.086 

-0.026 ± 0.060 

0.67 

0.985 

± 0.035 

0.006 ± 0.022 

0.76 

0.19 

0.625 

±0.181 

0.156 ±0.160 

0.41 

0.31 

Apodization 

1.009 ± 

0.088 

-0.001 ± 0.023 

0.98 

0.985 

± 0.026 

0.002 ± 0.008 

0.85 

0.58 

0.753 

± 0.097 

0.030 ± 0.042 

0.50 

0.91 

Bias-hardened estimators 
$ TT - M 11 

1.117 ± 

0.120 

-0.073 ± 0.054 

0.17 

0.968 

± 0.039 

0.022 ± 0.021 

0.30 

0.33 

0.451 

±0.136 

0.273 ± 0.098 

0.01 

0.11 

<j> TT - N 11 

1.068 ± 

0.116 

-0.024 ± 0.036 

0.52 

0.986 

± 0.036 

0.004 ± 0.014 

0.79 

0.79 

0.651 

±0.103 

0.074 ± 0.035 

0.04 

0.09 

fr TT - S 11 

1.082 ± 

0.118 

-0.038 ± 0.045 

0.43 

0.982 

± 0.038 

0.008 ± 0.017 

0.64 

0.62 

0.617 

±0.107 

0.107 ±0.045 

0.01 

0.06 


t Comparison is against the SMICA polarization-only reconstruction, rather than the MV reconstruction. 

^Comparison is against the SMICA TT reconstruction, rather than the MV reconstruction. 

This test only involves a change in the simulations used to perform the debiasing corrections and estimate error bars. We have quoted the scatter between the two sets of simulations, although this 
is an upper limit on the difference that we should actually see for the data. 









Table 3 Null tests for a set of lens reconstructions based on maps that should contain no CMB signal, or estimators that should be orthogonal to the lensing signal. The measured 
overall lensing amplitude A is given for the conservative (40<L<400) multipole range, as well as lower and higher multipoles. All should be consistent with zero. The PTE 
values are estimated by comparison of the y 2 values to the distributions from simulations. PTEs are given for the difference of the overall amplitude from zero, as well as for the 
X 2 of the difference from zero of the binned spectra. 


Null test summary values 


Low-L Conservative High-L 





VI 

oo 

< 

40 




40 

< 

L< 400 




400 

< 

L < 2048 




A 

£ 

+1 



PTE i=0 


A 

± O'A 



PTE ; i=0 

PTE bins 


A 

£ 

+1 



PTE i=0 

PTE bins 

Half-ring difference 

(—0. f 16 

± 

0.371) 

X 

io - 2 

0.78 

(-0.331 

+ 

0.154) 

X 

10- 2 

0.03 

0.04 

(-0.214 

± 

0.434) 

X 

10- 2 

0.54 

0.18 

Half-mission difference 

(-0.282 

± 

0.565) 

X 

io - 2 

0.63 

(-0.215 

+ 

0.177) 

X 

10- 2 

0.21 

0.39 

( 0.075 

± 

0.545) 

X 

10- 2 

0.93 

0.83 


( 0.658 

± 

0.843) 

X 

io - 1 

0.43 

(-0.099 

+ 

0.290) 

X 

10- 1 

0.74 

0.46 

(-1.891 

± 

0.967) 

X 

10-' 

0.05 

0.41 


(-0.241 

± 

1.413) 

X 

io - 1 

0.87 

(-0.415 

+ 

0.428) 

X 

10- 1 

0.34 

0.27 

(-3.034 

± 

1.047) 

X 

10-' 

0.00 

0.03 

lf, EE 

( 1.138 

± 

1.075) 

X 

10° 

0.29 

(-0.266 

+ 

0.463) 

X 

10° 

0.59 

0.66 

( 7.748 

± 

13.27) 

X 

10° 

0.53 

0.56 

xjj TE 

( 0.214 

± 

0.463) 

X 

10° 

0.64 

(-0.249 

+ 

0.185) 

X 

10° 

0.18 

0.25 

(-3.183 

± 

1.938) 

X 

10° 

0.11 

0.75 

l/, EB 

( 0.484 

± 

0.655) 

X 

10° 

0.46 

( 0.057 

+ 

0.213) 

X 

10° 

0.79 

0.41 

( 1.014 

± 

1.222) 

X 

10° 

0.43 

0.97 

f TB 

( 1.352 

± 

2.563) 

X 

10° 

0.60 

( 0.378 

+ 

0.860) 

X 

10° 

0.63 

0.81 

( 6.597 

± 

5.144) 

X 

10° 

0.21 

0.24 
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5. Conclusions 

We have presented a first analysis of gravitational lensing in the 
Planck 2015 full-mission data set. Using temperature and polar¬ 
ization data, we make the most powerful measurement of CMB 
lensing to date, with an overall significance of around 40 cr. This 
paper is accompanied by the public release of the lensing poten¬ 
tial map (as well as simulations that can be used to characterize 
the map, e.g. for cross-correlation analyses with external data) 
and a conservative lensing bandpower likelihood. 

We have performed a number of tests to verify the inter¬ 
nal consistency of the lens reconstruction, in particular we have 
tested the compatibility between individual temperature and po¬ 
larization estimators, sensitivity to various analysis choices, and 
possible foreground contamination. 

Given Planck ’s unique full-sky coverage, the full-mission 
maps are unlikely to be displaced as the most powerful data 
set for cosmological lensing analysis for the next several years 
(until high signal-to-noise E- and fi-mode polarization measure¬ 
ments can be made on sufficiently large sky fractions). There 
is some room, however, for improvement on the analysis pre¬ 
sented here. The inverse-variance filtering we have used does 
not account for variations in the instrument noise level across 
the sky, which could be incorporated to reduce the reconstruc¬ 
tion noise (particularly in polarization, where instrumental noise 
dominates the signal). We have high-pass filtered all results to 
L > 8, due to concerns about the lensing mean-held correction. 
A more aggressive analysis could potentially remove this cut. 
The conservative L < 400 cut that we make in our baseline like¬ 
lihood, imposed because of mild evidence of a broad feature in 
curl-mode tests beyond this range, might also be relaxed in a 
future analysis, although only with minor improvements in the 
statistical power of the likelihood. 

The microwave background is a unique source for lensing 
studies, as the most distant and well-understood source plane 
that can be observed. The measurement presented here repre¬ 
sents an integrated measurement of the total matter distribution 
in the entire observable Universe, and another powerful test of 
the veracity of the concordance ACDM cosmological model. 
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Appendix A: Pipeline 

In this appendix we provide a more detailed technical overview 
of the pipeline that we use to extract lensing information from 
the Planck maps. Starting from input CMB sky maps the pipeline 
has the following three distinct analysis stages. 


(1) Filtering. In this step, we perform a linear operation on the 
maps, which performs three tasks. 

a. Mask out regions of the sky that may be contaminated by 
point sources or residual foreground emission (primarily 
from the Galaxy). 

b. Deconvolve the instrumental beam and pixel transfer 
functions. 

c. Downweight noisy modes to allow later steps of the 
pipeline to produce optimal (in the minimum-variance 
sense) estimates of the lensing potential and its power 
spectrum. 

Implementation details of the filtering procedure are de¬ 
scribed in Appendix A.l. The filtering step provides us with 
filtered CMB temperature and polarization multipoles de¬ 
noted by T (m , E lm , and B [m . 

(2) Quadratic Estimators. The filtered CMB multipoles are then 
fed into quadratic estimators designed to extract the statis¬ 
tical anisotropy that is induced by lensing. These estima¬ 
tors are formed by summing over pairs of CMB fluctuations 
with a weight function optimized to detect the off-diagonal 
contributions to their covariance matrix from (fixed) lenses. 
This step produces maps of the estimated lensing potential 
(j), and is described in detail in Appendix A.2. 

(3) Power Spectrum Estimation. To estimate the power spec¬ 
trum of the lensing potential, we exploit the non-Gaussianity 
induced by lensing. At first order in the power spectrum , 
lensing generates a contribution to the connected 4-point 
function (or trispectrum) of the observed CMB, which we 
probe using auto- and cross-spectra of the 2-point quadratic 
estimators 4> obtained in step (2) of the pipeline. There are 
two complications here that are worth noting: the auto- and 
cross-spectra of the lensing estimates have noise biases even 
in the absence of lensing, which must be carefully estimated 
and subtracted; and non-Gaussian contamination from unre¬ 
solved point sources can mimic lensing and bias the recon¬ 
structed power spectrum. Our approaches to these issues are 
described in Appendix A.3. 

We now proceed to describe the individual steps of the lensing 
pipeline in more detail. 

A.l. Filtering 

At this stage in the pipeline, we have noisy maps of the CMB 
for both data and simulations. To construct minimum-variance 
estimates of the lensing potential and its power spectrum, we 
need to filter these maps to downweight noise-dominated modes, 
as well as to deconvolve the transfer functions due to the finite 
beam size and pixelization. 

To derive the optimal (minimum-variance) filter to apply, we 
use a simplified model of the data given by 

di = ^ 3/ikSk + n h (A.l) 

k 

where di is a single long vector containing the pixelized T, Q, 
and U sky maps. The ‘pointing matrix’ J/yt takes an input sky 
signal (consisting of T, E. and B multipoles collectively indexed 
by k), convolves it with a beam and pixel transfer function, and 
then performs a harmonic transform to map space. The harmonic 
modes of the sky itself are denoted as s, t, and n, is the map noise 
realization. We have used indices i and k for quantities in map- 
space and harmonic-space respectively. Explicitly, the indices k 
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and i are mapped as 

k (le {T,E,B}J,m) (A.2) 

/ h» (Me{T,Q,U},p), (A3) 

where p denotes a pixel. In this notation, we model the pointing 
matrix as 

y <M,p),(x/m) = HfBfR^ f:m (h p ), (A.4) 

where H( is the HEALPix (Gorski et al. 2005) N s a e = 2048 pixel 

transfer function, B/ is the effective beam transfer function for 
the map under consideration, and R^ fm (h p ) represents the real¬ 
valued spherical harmonic for field X evaluated at the center of 
pixel p in map M. 

If both the sky signal and noise are Gaussian random fields 
(with covariance matrices (SkS* k ,) = and (<) = N, re¬ 
spectively) then the inverse-variance filtered estimates of the 
beam-deconvolved sky signal s are given in matrix notation by 
(Bunn et al. 1994) 

s = S“ 1 [S“ 1 +2/ r N- 1 J/]~ 1 (A.5) 


We use this equation to obtain the inverse-variance filtered T, 
E, and B multipoles that are fed into our quadratic estimators, 
with one additional rescaling discussed below. The large brack¬ 
eted matrix inverse is performed using conjugate descent with a 
multi-grid preconditioner, following Smith et al. (2007). 

For the purposes of filtering, we use a simplified model for 
the covariance matrices of su and in which they are diagonal 
in harmonic space and pixel space, respectively, with 

S — J:'m'~ ^XX’^lt'^mm'Q t , 

N = ( n(M,p)n(M',p ’)) = dMM'dpp'Np . (A.6) 

Here, Cf x are the fiducial theoretical CMB power spectra and 
Np gives the noise variance for pixel p in map M. With this 
diagonal assumption, we have ignored the cross-correlation C T ( E 
for the purposes of filtering. This is slightly sub-optimal, but has 
the advantage of allowing the temperature and polarization maps 
to be filtered independently. It also has the benefit of making 
analytical calculations for the estimator normalization simpler 
(see Sect. A.2). Also for simplicity, we take N' p to be constant 
over all unmasked pixels and given by 


N m = 

p 


10 800 


Apix 1 

4n M n 


(A.7) 


Here N™ is the map noise level in pK arcmin, A p ; x is the number 
of pixels in an Wide = 2048 HEALPix map, and M p is a mask 
map (which is zero for masked pixels, and unity otherwise). We 
always take the noise level to be the same for paired QJU maps. 
The mask term in Eq. (A. 7) takes the noise level to infinity for 
masked pixels. In the filtering operation, which only involves the 
inverse of the noise matrix, this sets masked pixels to zero. 

After evaluating Eq. (A. 5), we have a set of inverse-variance 
filtered sky multipoles sx.tm for X e {T, E, B). The simplified 
CMB+white-noise model above results in a slightly sub-optimal 
filtering procedure since it ignores the following: 


(1) variations in the noise level across the sky due to the un¬ 
even hit distribution (the Ecliptic poles receive more cover¬ 
age than the equator for the Planck scan strategy); 

(2) scale-dependence of the noise power; and 

(3) foreground power. 


In Planck Collaboration XVII (2014) we showed that the sub¬ 
optimality due to (1) for temperature-based lensing reconstruc¬ 
tion was small (accounting for it properly was estimated to im¬ 
prove the uncertainty on the overall lensing amplitude by 4%). 
The loss will be slightly smaller here for temperature, owing 
to a more even hit distribution in the full-mission data than in 
the nominal-mission data used in Planck Collaboration XVII 
(2014). The loss will be larger for polarization, where noise is a 
more important part of the total error budget, although we have 
not evaluated the size of the potential degradation. To compen¬ 
sate for (2) and (3) we make a post-correction by rescaling the 
modes by a quality factor Q x as a function of scale to obtain fi¬ 
nal filtered multipole that are closer to optimal. We denote these 
as X( m , and they are given by 


Xfm — Qp Sx. 


im* 


(A.8) 


These inverse-variance filtered multipoles are the input to the 
lensing estimators described in Sect. A.2. The choice of Q x is 
discussed below. 

For analytical calculations, it is useful to have an approxima¬ 
tion to the filtering procedure that is diagonal in harmonic space 
(i.e. neglecting the small amount of mode-mixing induced by 
the masking). In this approximation, the filtered multipoles are 
given by 

At m ~ FfX (m + h( m , (A.9) 

where n is a noise realization and 


F 


x 

t 


f^XX 


Qt 


+ N XX ' 


(A. 10) 


Here, C xx are the fiducial CMB auto-spectra used in Eq. (A.6), 
and the quantity 


N xx = (HcBtT 2 3 




10 800 J 


(A.11) 


is the pixel- and beam-deconvolved noise power spectrum for 
the field X, where N^ v is the appropriate map white-noise level 
in /rK arcmin. We choose Q x such that 


f-1 
J sky 


(2 L + 1) 


2 \X em \ 2 * Ff, 


(A. 12) 


where f s ^ y = Y, P M p /Npa. is the unmasked sky fraction and the 
approximation sign indicates a smoothing over multipoles. 


A.2. Quadratic estimators 

Ensemble averaging over a fixed realization of the CMB lens¬ 
ing potential, the CMB covariance matrix acquires off-diagonal 
elements given by 

A<X fim| Z, 2m2 ) = g(-l)' w (^ ^ _ L M ^W x X L d,.M, (A. 13) 

with the fields X fm , Z (m e {Tf m , E fmi B ( m ). The covariance re¬ 
sponse functions TV XZ for the possible field combinations can 
be found in Okamoto & Hu (2003). They are linear in the 
(lensed) CMB power spectra, and are only non-zero if + {2 + L 
is even for even-parity combinations, e.g. TT and TE, and 
f\ + (o + B is odd for the odd-parity combinations TB and 
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EB. Furthermore, the r W xz are real for the even-parity com¬ 
binations and imaginary for the odd-parity. We use the off- 
diagonal covariance to estimate the lensing potential; after the 
inverse-variance-filtered sky maps are generated, pairs of fields 
X( m , Z( m e \ f f m , E( m , B/m 1 are fed into quadratic estimators to 
estimate the lensing potential. A general quadratic estimator is 
given as a function of input fields X and Z as 


xlm[X, Z] 


(~1) M y / t\ h L \ 

2 Z—k 1 iti\ m 2 -M I 

X (A.14) 


where e L is a set of weight functions that define the estimator 
x, along with the input inverse-variance filtered multipoles X, Z 
for a particular map or pair of maps. The choice of T, E, or B 
for X and Z depends on the specific estimator x. In most of what 
follows we shall drop the [X,Z] unless it is necessary to avoid 
ambiguity. 

Optimal lensing estimators use a matched filter for the 
lensing-induced covariance of Eq. (A. 13). We denote these as 
W xz for the fields X and Z. The temperature-only optimal esti¬ 
mator, for example, uses W TT = r W 7T | fid , where the ‘fid’ sub¬ 
script indicates that the lensing weight function is evaluated for 
the fiducial cosmological model given at the end of Sect. 2. For 
the odd-parity combinations, the matched filter is minus the co- 
variance response function, e.g. W TB = - TV 15 ^. Neglecting 
the lensing of any primordial fi-mode signal (certainly a valid 
approximation for scales f « > 100 with Planck polarization 
sensitivity), generally there are eight possible lensing estima¬ 
tors based on the possible combinations of T, E, and B. (This 
reduces to five estimators for our baseline analysis that uses a 
foreground-cleaned map of the CMB in temperature and polar¬ 
ization, and for which there is no distinction between TE and 
ET , for example.) 

The quadratic estimators defined by Eq. (A. 14) will receive 
contributions from non-lensing sources such as masking, beam 
asymmetry, and inhomogeneity of the instrumental noise. We 
determine and correct for this bias by averaging the reconstruc¬ 
tion xim over Monte Carlo simulations that include these effects, 
thereby estimating a ‘mean-field’ that we then subtract. 

The response of a quadratic estimator to the covariance given 
in Eq. (A. 13) (averaged over CMB realizations with a single 
mode (f>LM of the lensing potential held fixed) is, in the diago¬ 
nal approximation of Eq. (A.9), 


(x LM ) = <Plm, (A. 15) 


where 


K = 


---y w h i F f F r ■ 

2(2 L + 1) ' , J *1*2 L lll2 L c 1 12 


(A. 16) 


As a baseline, all of the results in this paper are obtained using 
the fiducial model at the end of Sect. 2 to calculate the TV XZ 
that appear in the estimator response. However, for cosmologi¬ 
cal parameter sampling with the lensing likelihood the normal¬ 
ization is recalculated using the appropriate CMB power spectra 
as the parameter space is explored. This is discussed further in 
Appendix C. 

Putting the above together, we form estimates of the lensing 
potential as 

&LM = ~ X LM )• (A. 17) 


We can also sum the individual estimators into a combined 
minimum-variance estimator (MV) as 





(A. 18) 


where the sum is taken over the eight lensing estimators (TT, 
EE, TE, TB, EB, ET, BT, and BE). 

Finally, we note an important implementation detail. 
Naively, evaluating an estimator of the form in Eq. (A. 14) re¬ 
quires <9(f^ ax L^ ax ) operations, which would be very computa¬ 
tionally expensive at Planck resolution. Fortunately, all of the 
weight functions used in this paper can be re-written as a sum of 
separable terms with the following form 
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This separability leads to estimators that can be evaluated with 
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(A.20) 


The separability of the weight functions also allows the 
response functions to be calculated quickly analytically; 
using the position-space approach given in Dvorkin & 
Smith (2009) they can be evaluated in <9((^ ax ) opera¬ 
tions. A sample implementation of these weights can be 
found at http://github.com/dhanson/quicklens in the file quick- 
lens/qest/lens.py. 


A.3. Power spectrum estimation 

Ensemble averaging over realizations of both the CMB lensing 
potential and the primary temperature and polarization fluctua¬ 
tions, the CMB becomes non-Gaussian. At first order in the lens¬ 
ing potential power spectrum cff this non-Gaussianity is mani¬ 
fest as a connected 4-point function (Hu 2001) 


(X{ imi Zf 2m2 C{ 3fn4 T){! 4m4 ) c — 2 ^ 
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{1 £2 L 
m\ m 2 M 
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(A.21) 


As might be expected given the presence of the covariance re¬ 
sponse functions TV in this expression, estimators to extract Cf 1 ’ 
from the trispectrum can be written using the quadratic estima¬ 
tors above as building blocks. We write the cross-spectrum of 
two quadratic estimators explicitly as 


L,xy 
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Fig. A. 1 Lens reconstruction bias terms for the MV lensing es¬ 
timator applied to the Planck 2015 SMICA maps. Dashed lines 
indicate regions where the bias term is negative. The fiducial 
ACDM theory power spectrum C^' fid is plotted as the black 
solid line. 


where y = £ p M p /N p \ x is the unmasked sky fraction, with 
M p the mask map used in Eq. (A. 7). We form estimates of the 
lensing potential power spectrum (based on estimators x and y) 
as 


I 

L,xy L,xy T-HnO 

- SCf I - A Cf I - A df I . 

L.^|ni Pmyic ^Ips 


(A.23) 


There are several correction terms here, which are discussed in 
more detail below. We also plot them for the MV reconstruction 
based on the Planck 2015 SMICA maps in Fig. A.l. 

The AC^Jno term represents the disconnected contribution 
to the 4-point function, which would be non-zero even in the ab¬ 
sence of lensing. We estimate this term by replacing some of the 
data fields with those from two sets of independent simulations 
(labelled MCI and MC2), and then average over realizations as 


ACfJ N0 = | - C^[A r MCl>Z M C2,CMC2jC ) MCl] 

+ C^ xy [X MC i,Z, Cmci>^] + C^ xy [X M c\,2,C,Dy lC i\ 

+ eg, [X , Z M c i , Cmc i . D] + 2mc i . C , D M c i ] 

- C*y [X M c i,Zmc2 ,CmcuD MC 2]\ . (A.24) 

/MCI, MC2 

This method for determining the disconnected bias can be de¬ 
rived from the Edgeworth expansion of the lensed CMB (Regan 
et al. 2010; neglecting the T(f> and E(f> correlations). If the sim¬ 
ulations have a covariance differing from reality by a linear (in 
general, anisotropic) correction, this form of the NO estimator 
using a mixture of data and simulations is insensitive to the cor¬ 
rection at linear order (Namikawa et al. 2013). The largest differ¬ 
ence we see between simulations and data is about 5% in power 
before correction, so the quadratic error should be at most about 
(0.05) 2 = 0.0025, corresponding to an error of less than 0.4 cr for 
aggressive bins at L < 400 rising to around 1.3 cr at L = 2000. 
Since we also match the data calibration, and add isotropic fore¬ 
ground and noise power to the simulations to match the observed 
amplitude and shape of the power spectrum of the data, the 


isotropic part of the error is actually close to zero. A residual er¬ 
ror only arises at second order due to our inability to capture any 
statistical anisotropy in the missing foreground and noise power, 
and hence should be much smaller than the above estimates. This 
makes the NO subtraction rather robust, despite it being substan¬ 
tially larger than the signal on small scales. Furthermore, the 
mixed form reduces correlations of the lensing power spectrum 
estimates between different multipoles (Hanson et al. 201 1) and 
with the measured CMB power spectra (Schmittfull et al. 2013). 

The AC^Ini term corrects for the non-Gaussian secondary 
contractions (the other permutations in Eq. A.21) due to lensing, 
as discussed in Kesden et al. (2003). We evaluate it using the 
flat-sky approximation in 2D Fourier space as 
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where t\+t 2 = (\ +t\ = L. Here the weight functions are flat-sky 
lensing weight functions; they can be found in Hu & Okamoto 
( 2002 ). 

We characterize any differences between the average lens¬ 
ing power obtained on simulations with the input power. This is 
determined by 
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(A.26) 


The MC correction is included to account for possible issues 
in the estimation of the normalization, mixing of power between 
power spectrum bins, and errors in the calculation of the N1 bias. 
Although in principle several of these effects should be included 
as multiplicative rather than additive corrections, we find that the 
MC correction is generally small enough (less than 10% of the 
fiducial lensing power spectrum) that it does not matter at our 
experimental sensitivity whether we apply the correction addi- 
tively or multiplicatively. 

Finally, we make a correction AC^Jps for the non- 
Gaussianity due to the shot-noise of unresolved point sources 
in temperature. This is given by 
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_ q 4 L, L 




(A.21) 


where S 4 is an estimate of the shot-noise amplitude and S 2 de¬ 
notes the weight function (non-zero only for TT) defined by 


W Ll 


(2 G + 1)(2G + 1)(2L + 1) 


ti { 2 L\ 

0 0 of- 


(A.28) 
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The R xS is given by Eq. (A. 16) but with c x t 2 L replaced by 
W( t f L - Following Osborne et al. (2014), the estimate for the 
shot-noise trispectrum amplitude is obtained in an analogous 
way to Eq. (A.23), but using the quadratic estimator formed with 
the point-source weight function above rather than the lensing 
weight function. The overall source amplitude is then estimated 
as 


S 4 = 


ZdiL+Dftp 2 ) 2 C s l 2s2 
£ L (2L+l)(7?f s2 ) 2 


(A.29) 


where the sums are taken over 100 < L < 2048. Note that R s2sl 
is given by Eq. (A. 16) but with both 'We^L and W} ( L replaced 

b y K 2 l- 

We use Eq. (A.23) to estimate the lensing potential power 
spectrum, which we then bin into bandpowers for plotting as 
well as constructing a lensing likelihood. We use Monte-Carlo 
simulations of the power spectrum estimation to construct er¬ 
ror bars for these bandpowers. There is one technical point here, 
which is that the realization-dependent AC^.Jno correction term 
can be cumbersome to compute. Instead, for the purpose of de¬ 
termining error bars, we use a semi-analytical approximation to 
this term given by 


AC' 


L,xy 


I _ 1 1 

\X Z C D\ — _ 

I NO, analytic L ’ ’ ’ 4(2 L +1) 


x X KUKe.ZfCf + (-l) A+4+i ^ iL Cf Cf), 

(A.30) 




where the empirical power spectra of the filtered fields are 


r xz 

u e 


J sky 

2L + 


j X X{mZ (m ■ 


(A.31) 


Appendix B: Bias-hardened estimators 

The quadratic lensing estimators of Appendix A.2 can be biased 
by non-lensing sources of statistical anisotropy such as beam 
asymmetry, noise inhomogeneity, and masking. In our fiducial 
analyses, we estimate these mean-field biases using Monte-Carlo 
simulations. However, this approach can be sensitive to the fi¬ 
delity of the simulations, particularly for the reconstruction on 
large angular scales where the mean-field biases can be orders 
of magnitude above the lensing signals of interest. As a cross¬ 
check, we also take a more data-dependent approach to cor¬ 
recting these biases. As with lensing, most sources of statistical 
anisotropy in a map can be associated with a covariance response 
function TVJ ^, and a held zlm that describes the spatial de¬ 
pendence of the anisotropy. It is possible to construct quadratic 
estimators for zlm , and then calculate and subtract the result¬ 
ing lensing bias. This procedure amounts to constructing spe¬ 
cial ‘bias-hardened’ estimators (Namikawa et al. 2013). Given 
a source of bias z, we construct a bias-hardened estimator with 
weight function 


<n4>~ 

W^~ Z - W l/ ’_— 

vv (i t 2 L ~ y hhL 


_ w~ 

(J^zz IxtjL ’ 


(B.l) 


where R? z and Rf are generalizations of Eq. (A. 16), and { L 
is the matched filter for TV^^. This method is potentially less 


susceptible to errors in the mean-held estimate (due to inaccu¬ 
racies in the simulations used to obtain it), although it has sev¬ 
eral caveats that are discussed in more detail in Namikawa et al. 
(2013). 


Appendix C: Power spectrum likelihood 


Here, we provide more detail on the components of our lensing 
likelihood (given in Eq. 3, and copied below for convenience) 

- 2log= Sf(Cf - Cf’*) [S -1 ]' 7 S L -{C^ - C£?' th ). (C.l) 
The bandpower binning function S E is given by 

£•00 > fidy_i 

Sf = -:----, (C.2) 

^^max fid^ y_l 

— ^min 

where Vl is an approximation to the variance of the power spec¬ 
trum estimate. For a given cross-spectrum between estimators x 
and y we use 

, = 2L+1 ^ 0 

2/ s - k ; * L 

which produces minimum-variance estimates of the lensing am¬ 
plitude for each bin in the limit that the sample variance of 
the lenses is negligible (a reasonable approximation for Planck, 
which is noise dominated on any individual mode). 

The dependence of the likelihood on cosmological parame¬ 
ters enters through the ‘theory’ spectrum C^' th . This is the ex¬ 
pected value of the estimated spectrum at each point in parame¬ 
ter space. It depends on the cosmological parameters 6 in several 
ways: 

• directly, through the theory spectrum C^L; 

• indirectly (but linearly) on C^\ g , via the theory-dependence 
of N m (Eq. A. 25); and 

• indirectly, and non-linearly, on the CMB power spectra 
Cj T , Cj E , and Cf E through the estimator normalization and 
N <1> (from the theory-dependent covariance response TV of 
Eq. A. 13). 


We neglect any other theory dependence, for example in the MC 
correction that is hard to quantify. In detail, for a given set of 
cosmological parameters 6, we should calculate 
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00,th 


CR^R^) I 

r L \0 --00 /r 

(Kfflf)| L ' e ^Ini, fid 


where the TV component of the response functions R in the nu¬ 
merator, as well as the N (l) term, are calculated using the CMB 
power spectra for parameters 6. 

Fully recalculating everything for each point in parameter 
space is prohibitively slow, but for small deviations from the 
fiducial model we can use a linearized approximation. The de¬ 
pendence on is already linear, and expanding the CMB 

power spectrum dependence to linear order about the fiducial 
model we can write 


C 


,00,th 


C' 


d MRf/Rtf) 


L \0 


dCp 


(TL-c;.|J 


c 


L Ifid 


I ^ 001 _ f00| \ 

1V1 LL',xy\^L’\f) ^L' Ifid) 


dAC 


L,xy 


N1 


dCj, 


(c {,I -d,I ), 

\ 1 le 1 Ifid / 

(C.5) 


25 










Planck Collaboration: Gravitational lensing by large-scale structures with Planck 


where j sums over the various CMB power spectra. The matrix 
(from Eq. A. 25) gives the exact linear dependence of N (1) 
on the lensing potential for fixed CMB power spectra, and can 
be pre-computed along with the other derivative matrices for the 
fiducial model. The binned likelihood of Eq. (C.l) then depends 
on a binned ‘theory’ power spectrum given in the linear approx¬ 
imation by 


S L C H>^ a S L 


■ Mf (c; 


i'\e 


c\ 


f I fid. 


)• 


(C. 6 ) 


where now a sums over both the C® and CMB power spectra 
terms, and M“' c can be pre-computed in the fiducial model. 

The Planck lensing likelihood consists of bandpower esti¬ 
mates and the covariance Z, ; , a set of binning functions 

Sf 1 , a set of linear-correction kernels M" J , and the band values 
of the linear-correction kernel applied to the fiducial spectrum 
Mf l C“\ m . Note that since the lensing estimator is independent 
of the CMB calibration, for a theory model C J ( and map calibra¬ 
tion parameter y ca i, the theory CMB power spectra that are used 
in Eq. (C. 6 ) are those matched to the data, i.e. C ] ( \g = C J c /y ^ al . 

When evaluating ‘lensing-only’ likelihoods, where the CMB 
power spectra are only very weakly indirectly constrained, we 
fix the Cj\o used in the linear correction to a best-fit to the full 
data. When calculating joint lensing and CMB results, the cor¬ 
rection depends on the CMB power spectra at each point in 
model space. 


Appendix D: CIB model 

The cosmic infrared background (CIB) is a diffuse sky sig¬ 
nal which begins to dominate over the CMB at frequencies 
v > 300 GHz. It is generated by dust that is heated by UV light 
from young stars and then re-radiates thermally in the infrared. 
The CIB contains approximately half of the total extragalactic 
stellar flux, and has excellent redshift overlap with the CMB 
lensing potential (Song et al. 2003). In Sect. 3.2 we use the 
CIB as a tracer of the lensing potential to estimate the lensing- 
induced B-mode signal. This requires a model for the cross¬ 
correlation between the CIB and the lensing potential. We cal¬ 
culate this in the Limber approximation as (Limber 1954) 

C 0 = f d X - 2 K™(x)Kl(x)P{k = L/at; m ~x), (D-D 

a X 


where P(k\ r/) is the (non-linear) matter power spectrum at co¬ 
moving wavenumber k and conformal time rj. The lensing kernel 
in flat models is given by 


<(x) = ~ 


3 Q m H 0 2 
a(x) 



(D.2) 


where a(x) is the scale factor at conformal time tjo — X and X * 
denotes the last scattering surface at z ~ 1100. To model the 
CIB, we use the simple SSED model of Hall et al. (2010), which 
has a kernel given by 


= (D ' 3 > 

where b c is an overall normalization and z c - cr z = 2 describe 
the redshift distribution of the CIB intensity. The redshift z is 
evaluated at conformal time 770 - X - The spectral energy distri¬ 
bution (SED) of a typical CIB source is described by /„. It is 


modelled as a modified blackbody with temperature T = 34 K, 
and spectral index [5 = 2 up to a threshold v’ ~ 4955 GHz, where 
it transitions to a power-law decay with index a — 2: 


|[exp(^) - i] V +3 (v < v0, 
"* V ” \[ex P (|f) ~ l] * v '^ +3 ( 7 ) “ (v>A 


(D.4) 


The precise value of V is chosen to match smoothly the slope 
of f v on both sides of the transition frequency. In this model, the 
modified blackbody component of the SED is generally what 
we observe; even at an observation frequency of 1200 GHz (i.e. 
250 pm), the power-law transition does not occur until z ~ 4. Our 
theory curve is plotted in Fig. D.l, as well as the measured cross¬ 
correlation of the Planck 2015 545 GHz map with the MV lens¬ 
ing potential estimate (over 40 % of the sky). In addition we have 
plotted the bandpowers from the rigorous analysis presented in 
the Planck 2013 data release (Planck Collaboration XVIII2014). 
We have set the normalization b, for the theory curve such that 
the overall amplitude is consistent with both measurements; cal¬ 
culating x 2 values against the theory curve we find PTEs of 73 % 
for the MV lensing potential estimator and 49 % for the for the 
13 plotted bins we find a X 2 PTE with respect to the theory curve 
of 49%. 


S-( 
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s 
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Fig. D. 1 Cross-correlation between the MV lensing potential es¬ 
timate and the Planck 545 GHz channel (red points). The cross¬ 
correlation uses a Galaxy mask leaving 40 % of the sky. The ear¬ 
lier result from Planck Collaboration XVIII (2014) using a TT 
lensing estimator is shown with green points. The simple CIB 
model discussed in Appendix D is plotted as the black dashed 
line. 


Appendix E: Parameter dependence of in 
ACDM models 

In this appendix we discuss the dependence of the lensing poten¬ 
tial power spectrum at multipoles L > 100 on the parameters of 
ACDM models, providing results that are used in the discussion 
of the ‘lensing-only’ constraints in Sect. 3.5. Since these con¬ 
straints adopt a prior on n s that is rather tighter than the resulting 
posterior widths for D m /r and /;, we ignore dependencies on n s 
throughout this appendix. We caution the reader that for joint 
constraints, or fits to the CMB power spectra alone, fractional 
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variations in n s can be comparable to those in Q m /z 2 and h, so 
the dependence of Cf P on n s cannot be ignored. The dependence 
of Cf on the matter density was recently discussed in Pan et al. 
(2014); we follow their discussion closely. 

For multipoles L > 20, the angular power spectrum of the 
lensing potential can be Limber-approximated as 

L 4 cf = 4 r dx (k 4 Pv)(k = L/ X \ TJ 0 -x)0 -x/x.f (E.l) 
Jo 

in spatially-flat models. Here, Py(k; 77 ) is the equal-time dimen¬ 
sional power spectrum of the Weyl potential M J . Well after baryon 
decoupling, we have 1 

k 4 Pv(k- 77 ) oc k<P K (k)T\k)g 2 (a) , (E.2) 


where ViJk) is the dimensionless power spectrum of the primor¬ 
dial curvature perturbation R, T(k) is the transfer function, and 
g(a) is the growth function for the Weyl potential at scale factor 
a normalized to unity at high redshift. In ACDM, 


50 m E(a) r da' 

2 a Jo [a'E(a')] 3 ’ 


(E.3) 


where Q m is the matter fraction at the present and E(a) = 
H(a)/H d . For O m = 0.3, the current value of the growth function 
is 0.78 and t/lng/t/lnQ m = 0.23. Close to O m = 0.3, we can 
approximate the cosmology dependence of the growth function 
at redshift zero as g(a = 1) « 0.78(Q m /0.3) a23 . For L > 100, the 
lensing potential mostly probes early enough times that g(a ) = 1 
is a good approximation (Pan et al. 2014). However, if we wish 
to express the amplitude of P^(k) at early times in terms of cr 8 
(at the present), we require g(a = 1). 

The transfer function is normalized to unity for k <sc k eq , 
where k eq ~ 10 2 Mpc the inverse of the comoving Hubble 
parameter at matter-radiation equality. For k » k eq we have 
T(k) ~ (k/k eq y 2 \n(k/k eq ). Keeping the radiation energy den¬ 
sity fixed today, k eq oc cj m , where w m = O m h 2 . As well as k eq , 
the other scales that enter the transfer function are the sound 
horizon at the drag epoch (when baryons dynamically decou¬ 
ple) r d rag , and the diffusion scale k D . The broad-band shape is 
controlled by k eq , while rd rag determines the scale below which 
baryon suppression effects are important (see e.g. Eisenstein & 
Hu 1998). The non-zero baryon fraction Wb/« m tends to sup¬ 
press the transfer function on small scales since (essentially) un¬ 
clustered baryons impede the growth of the CDM density pertur¬ 
bation up to the time after decoupling when baryon infall into the 
CDM potential wells is complete. The further effect of baryons 
is to imprint baryon acoustic oscillations in the transfer function 
for kr dl ag > n and k < k V) . Baryon acoustic oscillations have 

little impact on since they are washed out by the line-of- 
sight integration. We therefore ignore them here and, where an 
explicit form for the transfer function is required, for illustrative 
purposes we use the ‘no-wiggle’ fit of Eisenstein & Hu (1998). 
For the range of parameters of interest, we can assume that the 
physical baryon density (±> d is fixed, so the transfer function can 
be thought of as a function of k/k eq , kr dmg and u m . 

Ignoring departures of the growth function from unity and 
approximating V'/Jk) as scale-invariant with amplitude A s , we 
can write 


L 4 C 7 oc A,4q 


/' 


dx 


( eq-T 


L 

X X* 


(1 


(E.4) 


where x = xix* an d 4q = k eq x* ~ 140. The A s 4q prefactor 
has a simple physical interpretation: CMB photons are typically 
deflected by { eq lenses of size k~ q , and each mean-squared de¬ 
flection is proportional to A s , giving a mean-squared deflection 
proportional to A s ^ eq . Given the ^-dependencies of the trans¬ 
fer function through k/k eq and kr^, and its direct parameter- 
dependence on (i) m , we see that in ACDM, the lensing power 
spectrum at L > 100 is (mostly) determined by only A s , L/t cq , 
L/T’drag, and a> m . Here, the multipole £di- ag = X* / r dra g character¬ 
izes the BAO scale at last-scattering. The integral in Eq. (E.4) 
determines the shape of the lensing power spectrum. Defining 



^drag 


=/' 


7 — ,o> m = dx — \T 2 — Ul-xj 


k& q x 


X* x 


(E.5) 


the dependence on the cosmological parameters (e.g. ( eq ) in the 
vicinity of the fiducial model is, at any L, approximately a prod¬ 
uct of power laws with exponents 


7 z^ q = - d In fiJd ln(L/4 q )| nd , 

(E. 6 ) 

T ag = - dlnf L /dln(L/e dlag )\ M , 

(E.7) 

zz™ = dln fi/d Inco m | nd , 

(E. 8 ) 


each evaluated for the fiducial cosmology. Note also that n and 
m 4k ' s determine the local power-law slope of C^, with 


n \; q + n d ™ g = -d\n{L A C^)/d\nL. (E.9) 


On the angular scales where Eq. (E.4) is valid, only modes with 
k > k eq contribute. For these scales, ( k/k eq )T 2 (k ) increases with 
k eq giving tf 4 > 0 (and monotonic increase of if 4 with L). 
Around L = 200, the scales that dominate the S/N of a lens¬ 
ing amplitude measurement from Planck, ~ 0.9. This, and 
A s , dominate the parameter dependence of C,f''. The effect of 
baryon suppression is reduced as w m increases at fixed Wb, so 
77™ > 0 and also rises monotonically with L. Around L = 200, 
we find 77 ™ ~ 0.3, which gives a small, but non-negligible, fur¬ 
ther parameter dependence to Cf P . Finally, the dependence on 
^dr ag is very small, since n'f dg < 0.15 (with peak value around 
L ~ 100, roughly equal to £dmg)- Putting these pieces together, 
we find a dependence of L 4 C‘x on cosmological parameters scal¬ 
ing approximately as 


l 4 c 00 



(E. 10) 


for L > 100 , where the dependence on 7 ( i rag is very weak. 

We can express the parameter dependence of C *' 4 in terms 
of the usual ACDM parameters. The comoving distance to last¬ 
scattering, x*- is 


X* 



da 

a 2 H(a) ’ 


(E. 11) 


where a„ is the scale factor at last-scattering (whose dependence 
on cosmology can be ignored here). At fixed radiation energy 
density today, the dependence on Q m and h is approximately 


0 . 4,-1 


X* - EC 4 /* 


(E. 12) 


around Q m = 0.3 and li = 0.77 This differs from the (D. m h 2 )~ lp2 
of an Einstein-de Sitter universe due to the late-time effects of 
A. Using k eq oc Q m /z 2 , we have 


t We ignore the small effects of non-zero neutrino mass in the ( eq oc O^/z. (E. 13) 

base ACDM model, so that the linear growth of structure is scale- - 

independent well after baryon decoupling. * We ignore the small dynamical effects of radiation at high redshift. 
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At fixed baryon density, the sound horizon at the drag epoch 
scales like rd rag or oj “ 0 ' 25 in the vicinity of the fiducial model, so 


that 


P rv O - 0.15 0-5 

^drag ^ 

(E.14) 

Using these dependencies in Eq. (E.10), we find 


L 4 Cf cc A s (n 0 m 6 li) nL+1 o 0 m I5xn * a8 - a2x ' ! ” , 

(E.15) 

where 


til = - 0.5 x n 4lJS + 2n™ 

(E.16) 


is approximately 1.5 at L = 200. The remaining dependence on 
£2 m is very weak, with 0.15xn* dg -0.2xn™ » -0.05 at L = 200. 
Equation (E.15) gives the main dependence of the lensing power 
spectrum on the cosmological parameters of ACDM models at 
multipoles L > 100. 

E. 1. Normalization by cr 8 

The dependence of the lensing power spectrum on A s can be 
eliminated in favour of cr 8 . We use the Poisson equation to relate 
the dimensional matter power spectrum at a — 1 to that of the 
Weyl potential, giving 

P s (k- a= 1) oc u)^k A Pv(k\ a = 1). (E.17) 

It follows that the variance of the density contrast in spheres of 
radius R is 

cr 1 2 oc a> m 2 g 2 (a = 1 )J dkk 2 P K (k)T 2 (k)W 2 (kR) 

J r*oo 

dxx 3 T 2 (xk eq )W 2 (xk eq R) , (E.18) 
o 

where we have assumed that Pnik) is scale-invariant in passing 
to the second line. The window function 

W(kR) = 3/ ' {kR) = cos kR - sin kR) , (E. 19) 

kR (kR ) 3 

is the Fourier transform of a normalized spherical top-hat func¬ 
tion of radius R\ following convention, we take R = 8/z 1 Mpc. 

The integral in Eq. (E.18) is a function of k eq R, k cq r iirdgl and 
co m at fixed baryon density (with the latter two from the trans¬ 
fer function). Evaluating the logarithmic derivative with respect 
to these parameters at the fiducial values, using the ‘no-wiggle’ 
transfer function from Eisenstein & Hu (1998), we find that the 
integral scales as (k eq Ry lA u )\[; 45 , with only a very weak further 
dependence on A cq r c i rag , in the vicinity of the fiducial model. 
Finally, the parameter dependence of <x? is 

OC A s keq(U”V(a = IX^eq/T 1 ) -1 ’ 4 ^ 45 

oc AsQ'jf/i 3 - 5 . (E.20) 

This parameter dependence agrees well with direct finite- 
difFerencing of <x 8 calculations with CAMB. 
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